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FOREWORD 


The  resoziTOh  work  presented  in  this  report  was  performed  by  Bolt, 

Beranek  &  Hewman,  Inc.,  Cambridge,  Massachusetts,  for  the  Flight  Dynamics 
Laboratory,  Directorate  of  Aeromechanics,  Deputy  Consnander/Tschnology, 
Aeronautical  Systems  Division,  Wright-Patteraon  Air  Force  Base  under  AF 
Contract  fh:  AF33(6l6)-7973*  This  research  is  part  of  a  continuing  effort 
to  provide  design  devolopnent  information  and  prediction  techniques  for 
structural  vibrations  and  dynamics  in  flight  vehicles  which  is  part  of 
the  Air  Force  Systems  Command's  Applied  Research  Hrogram  75CA,  the  Mechanics 
of  Flight.  The  Project  Nr  is  I370,  "Dynamic  Problems  in  Flight  Vehicles" 
and  the  Task  Nr  is  137009, "Methods  of  Vibration  Prediction  and  Control*. 

Lt.  K.  A.  Wrobel  of  the  Flight  Dynamics  Laboratory  was  the  Project  Engineer. 
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effort  under  the  aboTW  mentioned  contract.  The  work  presented  in  this 
report  represents  a  joint  effort  of  the  authors,  with  Dr.  E,  E,  linger 
acting  as  project  supervisor. 

The  authors  gratefully  acknowledge  the  assistance  rendered  by 
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ABSTRACT 


New  methods  are  outlined  for  dealing  with  the  vibration  responses 
of  complex  flight  vehicle  structures  to  local  and  to  diffuse  acoustic 
excitation.  Energy  absorption  at  structural  Joints  and  acoustic  radi¬ 
ation  resistance  are  shown  to  be  important  In  establishing  levels  of 
these  reaponces.  Some  experimental  results  pertaining  to  energy  ab¬ 
sorption  coefficients  and  radiation  resistance  are  given,  and  proce¬ 
dures  for  estimating  the  latter  are  discussed. 

Feasibility  studies  of  vibration  absorbers  utilizing  viscoelastic 
spring  elements  and  distributed  mass  systems  and  of  vibration  isolators 
composed  of  viscoelastic  leaf  springs  are  sununarlzed.  Only  the  latter 
are  found  to  possess  some  practical  advantages  over  conventional  sys¬ 
tems. 

The  results  of  experiments  are  presented  which  demonstrate 
that  sound-to-structure  coupling  may  be  reduced  significantly  by  the 
use  of  beams  of  special  design  (See  page  28)  whose  stiffness  de¬ 
creases  with  increasing  frequency.  An  analytical  investigation  is 
summarized  which  shows  that  generally  damping  of  only  the  plates 
of  beam-plate  systems  may  be  more  desirable  than  damping  of  only 
the  beams. 


Results  of  exploratory  studies  are  presented  dealing  with  the 
mechanisms  responsible  for  the  damping  of  aircraft  structural  Joints, 
and  analytical  investigations  of  the  relation  between  structiiral  Joint 
absorption  coefficients  and  loss  factors  are  summarized.  Lo3s  factor 
data  pertaining  to  two  actual  aircraft  substructures  are  appended. 
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INTRODUCTION 


The  increasing  abundance  of  Intense  sources  of  broadband  vibra¬ 
tion  asaocjafced  with  present  and  future  air/space  veblcles  tends  to 
make  control  of  structural  vibratory  responses  essential  if  struc¬ 
tural  fatigue  failures  and  equipment  malfunctions  are  to  be  kept 
within  reasonable  limits. 

The  design  engineer  has  available  analytical  methods  and  tech¬ 
nological  developments  that  permit  him  to  design  simple  structures 
to  withstand  slmp^lc  (narrow-band)  excitation.  However  no  methods 
appear  available  at  present  that  permit  him  to  deal  similarly  with 
complicated  structures  subject  to  more  complex  excitation. 

The  work  summarised  in  this  report  was  undertaken  as  a  step 
toward  providing  the  lacking  design  information.  Part  I  of  this 
report  outlines  studies  directed  at  identifying  the  parameters  that 
affect  the  responses  of  complicated  structures  most  directly:  Sec¬ 
tion  lA  deals  with  systems  excited  at  a  point.  Section  IB  with  sys¬ 
tems  excited  by  a  diffuse  acoustic  field.  Part  II  deals  with  more 
specific  vibration  control  techniques.  Sections  IIA  and  IIC  describe 
feasibility  studies  of  viscoelastic  leaf  spring  Isolators  and  non- 
classical  vibration  absorbers.  Section  IIB  demonstrates  the  reduc¬ 
tion  of  sound  to  structure  coupling  attainable  with  beams  designed 
to  have  decreasing  flexural  rigidity  with  increasing  frequency.* 
Results  of  an  analytical  study  related  to  optimum  placement  of 
damping  treatments  are  discussed  in  Section  IID  and  results  of 
experiments  aimed  at  revealing  the  salient  features  of  the  damping 
inherent  in  aircraft  structural  Joints  are  presented  in  Section  HE. 
The  final  section  IIP  susnirarizes  Initial  analytical  results  per¬ 
taining  to  the  relation  between  structural  loss  factors  and  absorp¬ 
tion  coefficients.  Appendix  I  constitutes  a  detailed  discussion  of 
modal  densities:  Appendix  II  summarizes  loss  factors  measured  on 
two  substructures  of  an  actual  aircraft. 


*These  beams  are  the  subject  of  patent  applications  being  filed  in  the 
United  States  and  other  countries  by  Bolt,  Beranek  and  Newman.  Inc. 


Manuscript  released  by  the  authors  30  March  1962 

for  publication  as  an  ASD  Technical  Documentary  Report. 
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SECTION  I 


PARAMETERS  AFFECTING  RESPONSES  OF  COMPLICATED  STRUCTURES 


If  one  attempts  to  solve  in  detail  the  various  coupled  equa¬ 
tions  of  motion  that  describe  the  dynamic  behavior  of  a  structure 
one  encounters  Intractable  mathematics  for  all  but  the  simplest 
structures.  One  may  therefore  wish  to  sacrifice  detail  for  tract- 
ability  and  search  for  a  method  that  will  permit  one  to  describe 
In  simple  terms  the  average  vibration  of  a  complex  structure,  so 
that  the  salient  parameters  of  the  structural  response  may  be 
understood  and  controlled  In  design. 

A.  Point-Excited  Structures;  Room  Acoustics  Analogy 

In  describing  the  acoustics  of  a  room  one  faces  many  of  the 
problems  that  one  encounters  also  In  determining  the  vlbraulons 
of  complicated  elastic  structures.  In  room  acoustics  one  cannot 
solve  the  wave  equation  in  detail  because  of  the  complexities  of 
geometry  and  of  the  Interactions  between  the  waves  and  the  room 
boundaries.  Instead,  one  recognizes  the  diffuse  character  of  the 
sound  field  and  may  then  calculate  some  of  its  average  character¬ 
istics,  based  on  some  average  properties  of  the  room  and  its 
boundaries. 

Similarly,  detailed  solutions  of  the  equations  of  motion  of 
structural  systems  get  out  of  hand  very  rapidly  as  the  complexity 
of  these  systems  ^  ncx*Ccis  es.  Therefore  one  is  led  to  emulate  room 
acoustics,  postulate  a  diffuse  wave  field  in  the  structure  to  be 
studied,  and  attempt  to  calculate  mean  properties  of  the  vibration 
from  relatively  rough  descriptions  of  the  structure  and  its  boun¬ 
daries. 

The  work  reported  In  this  section  attempts  to  set  down  the 
reasoning  as  applied  to  complex  structures  (composed  of  beam  and 
plate-like  sub-structures),  to  delineate  the  important  parameters, 
and  to  describe  some  experimental  results  that  validate  the  theo¬ 
retically  developed  relations. 

1,  Basic  Concepts 

The  energy  E  with  which  a  structure  of  total  mass  M  vibrates 
may  be  expressed  as 

E  =  M  v^  (A.l) 

in  terms  of  the  mean  square  velocity  v  of  the  structure,  if  it  is 
assiimed  that  this  energy  is  uniformly  distributed. 
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If  there  are  mechanisms  that  remove  energy  from  the  structure 
(as  there  alwa3rs  are  In  practical  systems)  then  one  must  supply 
energy  to  the  structure  in  order  to  maintain  a  given  level  of  vi¬ 
bration.  When  the  power  supply  is  "ttarned  off”,  the  energy  in  the 
structure  may  be  expected  to  decay.  If  this  decay  is  assumed  to 
vary  as  e~°^,  where  t  denotes  time,  then  the  power  loss  is 
-dE/dt=6E.  In  the  steady  state  the  power  input  P  must  equal  the 
power  loss,  so  that  one  may  write 

p  =  5  Mv^  =  Mv^  .  (A.2) 

The  reverberation  time  T=13o8/b  Introduced  here  is  a  quantity  that 
lends  itself  well  to  experimental  determination;  it  is  defined 
(ref.  l)  as  the  time  required  for  the  energy  to  decrease  by  a 
factor  of  10°. 


At  frequency  f  the  power  loss  in  the  interior  of  the  struc¬ 
ture  may  be  described  (ref.  2)  in  terms  of  the  structural  loss 
factor  T),  as  P.  .=27rfTiE.  Losses  at  the  boundaries  of  the  struc¬ 
ture,  on  the  otnir  hand,  may  be  more  conveniently  treated  in  terms 
of  the  concept  of  absorption  coefficients.  The  absorption  coeffi¬ 
cient  of  a  boundary  la  defined  as  the  ratio  of  absorbed  power  to 
incident  power  at  that  boundary.  If  a  system  of  waves  moving 
along  a  structure  encounters  a  boxmdary,  some  energy  is  absorbed, 
the  rest  is  reflected  and  propagated  until  it  meets  another  boun¬ 
dary,  etc.  The  energy  absorbed  at  one  encounter  with  the  boundary 
is  yE,  where  7  denotes  the  absorption  coefficient  as  defined  pre¬ 
viously.  The  munber  of  boundary  encounters  per  unit  time  may  be 
expressed  as  where  e  denotes  the  group  velocity  of  the 

waves  on  the  structure*  (the  velocity  at  which  energy  propaigatea) 
and  L  the  mean  free  path  (the  mean  distance  travelled  by  the  wave 
traln*^etween  boundary  encounters).  Hie  mean  power  loss  per  unit 
time  at  the  boundaries  thus  Is  where  y  denotes  the 

mean  absorption  coefficient.  ®  “ 

If  the  boundary  is  not  xanlform,  y  may  be  obtained  from 

yL  =Zy^L^  (A.3) 

where  7-  is  the  absorption  coefficient  of  the  i—  piece  of  boun¬ 
dary  (which  has  length  L. ),  and  where  L=2  L.  denotes  the  total 
boundary  circumference.  ^By  combining  the  relations  for  Internal 
and  boundary  power  losses  one  may  obtain 


♦For  a  plate  of  flexural  rigidity  D  and  surface  density  p,. 
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(A.4) 


P  =  6E  =  [2TrfTi  +  YCg/\,]E 


Kosten  (ref.  3)  has  shown  that  for  two-dimensional  systems  the  mean 
free  path  I^=TrS/L,  where  S  is  the  total  surface  area.  Hence  the 
foregoing  result  may  be  written  as 

6  =  =  2TffTi  +  ^  L,  ,  (A.5) 

which  is  closely  related  to  a  result  given  by  Cremer  (ref,  4). 

Eqs.  (A. 2)  and  (A.5)  contain  the  salient  results  of  the  foregoing 
discussion.  Eq.  (A, 2)  permits  calculation  of  from  Input  power, 
mass  M,  and  reverberation  time  T  (or  decay  constant  6),  and  Eq. 

(a.5)  expresses  T  or  6  in  terms  of  Internal  losses  (described  by  t]) 
and  boundary  absorption  coefficients  y. 

The  assumptions  underlying  the  foregoing  derivations  are; 

1.  Energy  uniformly  distributed  over  structure  (wave  field 
diffuse;  wavelength  much  smaller  than  mean  free  path). 

2.  Exponential  decay  of  energy  with  time, 

3.  droup  velocity  (not  phase  velocity)  ind  Kosten’ a  mean 
free  path  relation  apply. 

4.  Boundary  absorption  contributions  add  linearly  (Eq,  A. 3). 

The  general  validity  of  these  assumptions  is  by  no  means 
self-evident.  The  results  of  some  experiments,  undertaken  to 
verify  applicability  of  these  assumptions,  are  discussed  in  the 
following  section. 

2,  Experimental  Verification 

1)  The  validity  of  Eq,  (A. 2)  and  of  the  underlying  assumption 
of  uniform  energy  distribution  was  verified  by  measuring  the  Input 
Impedance  of  a  point-driven  plate  and  comparing  the  measxired  values 
to  those  predicted  from  relations  involving  Eq,  (A, 2)  and  separately 
measured  reverberation  times.  Good  agreement  was  obtained  over  the 
entire  test  frequency  range. 

2)  Exponential  decay  of  energy  was  verified  for  a  number  of 
beam  and  plate  combinations  (excited  by  impacts)  by  means  of  a 
graphic  level  recorder.  (Dils  recorder  essentially  plots  log  E  vs 
time,  so  that  a  straight  line  record  corresponds  to  exponential 
decay.  Virtually  straight  lines  were  obtained  In  all  experiments 
of  this  type. 
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3)  Applicability  of  group  velocity  and  of  Kosten's  mean  free 
path  expression  was  verified  by  means  of  an  experiment  which  pro¬ 
vided  a  boundary  portion  of  unity  absorption  coefficient j  i.e.j  a 
means  for  permitting  all  incident  energy  to  leave  the  structure 
under  consideration.  A  plate  was  used  in  these  experiments  as  the 
primary  structiire,  a  second  highly  damped  plate  attached  to  the 
first  over  a  portion  of  one  of  its  edges  served  to  provide  a  boun¬ 
dary  with  unity  absorption.  This  "known"  absorption  coefficient 
was  compared  with  values  obtained  experimentally  and  computed  from 


(A.6) 


In  this  relation,  which  may  be  obtained  from  Eq,  (A. 5)*  T  denotes 
the  reverberation  time  obtained  with  7^=0  (l.e.,  before  adding  ab¬ 
sorbing  structures  at  the  i£il  portion  of  the  boundary),  and 
denotes  the  reverberation  time  observed  after  addition  of  the 
absorber  extending  over  Equation  (A, 6)  is,  in  fact,  useful 

for  the  practical  evaluation  of  absorption  coefficients  in  general. 


Again,  good  agreement  between  measured  7.  and  the  predicted 
value  of  unity  vfas  observed,  within  the  limitations  of  the  experi¬ 
ment. 


4)  The  linear  additive  property  of  absorption  coefficients 
indicated  in  Eq,  (A, 3)  was  verified  by  a  number  of  experiments 
involving  plates  with  boxandary  absorption  provided  by  various 
lengths  of  damping  tape  and  screwed-on  beams,  variously  oriented. 
The  additive  property  was  found  to  hold  for  all  configurations 
tested. 

3.  Results 


a)  Limits  of  Validity 

The  proposed  equations  should  not  be  applied  blindly  to  all 
situations.  One  may  readily  visualize,  for  example,  that  absorb¬ 
ing  devices  placed  near  nodal  positions  of  a  structure  will  have 
little  effect  on  the  motions  of  the  structure,  whereas  the  absorb 
tlon  of  devices  located  near  anti-nodal  positions  may  be  greatly 
enhanced. 

In  addition,  absorbing  structures  that  are  small  compared  to 
a  flexural  wavelength  on  the  primary  structure  may  have  their 
absorption  decreased  because  they  can  not  interact  effedtlvely 
v/ith  the  wave,  and/or  they  may  have  their  absorption  Increased 
somewhat  due  to  effects  analogous  to  acoustic  diffraction  (refs, 
5,6). 
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b)  Absorption  Coefficients* 

i)  Summary  of  Experimental  Results 

Figures  1  and  2  summarize  the  results  of  absorption  coeffi¬ 
cient  measurements  on  aluminum  plates  of  l/l6  and  1/32  inch 
thickness,  respectively.  The  various  absorbers  that  were  used 
are  also  shown  schematically  in  these  figures. 

In  analogy  to  room  acoustics  one  may  assign  structural 
absorbers  to  three  categories; 

1)  High  frequency  absorbers,  whose  absorption  coefficient 
Increases  with  frequency,  e,g.,  damping  tapes. 

2)  Mid-frequency  absorbers,  whose  absorption  peak  occurs 
somewhere  in  the  mid-frequency  range  (say,  between  300  and  8000 
cps)  and  is  governed  by  a  resonance  phenomenon,  e.g.,  bolted-on 
beams. 

3)  Low  frequency  absorbers,  whose  absorption  is  highest  at 
low  frequencies  and  decreases  to  very  small  values  at  high  fre¬ 
quencies,  e.g,,  parallel  beams  with  dissipative  material  on  the 
plate  between  them. 

ii)  Absorption  Mechanisms 

The  previously  stmanarlzed  experiments  were  undertaken  pri¬ 
marily  to  demonstrate  the  absorption  coefficient  concepts.  The 
mechanisms  that  contribute  to  energy  absorption  need  further 
exploration,  but  some  conclusions  may  be  deduced  from  data  that 
are  available  so  far.  These  conclusions  are  as  follows; 

Damping  tapes  -  The  dissipative  action  of  these  tapes  is 
fairly  well  understood  (ref.  2;,  although  additional  work  Is 
required  to  interpret  available  analyses  in  terms  of  the  desired 
absorption  coefficients.  Theory  predicts  that  the  damping  action 
of  a  tape  is,  as  a  first  approximation,  proportional  to  the  ratio 
of  tape  foil  thickness  to  plate  thickness.  Absorption  coefficients 
measured  with  the  same  tape  on  plates  of  1/16  and  1/32  inch  thick¬ 
ness  were  found  tc  be  in  the  ratio  1:2  over  essentially  the  entire 
frequency  rangej  in  good  agreement  with  the  theoretical  prediction. 

Beams  bolted  to  plate,  with  viscous  Interlayer.  -  Energy 
losses  here  are  undoubtedly  due  to  deformation  of  the  interlayer 
and  thus  are  associated  with  relative  motion  between  the  beam  and 
plate.  Absorption  peaks  may  be  expected  at  maxima  of  this  rela¬ 
tive  motion,  l.e,,  at  resonance  of  the  plate  portions  between 
screws  (since  the  beam  probably  deforms  and  moves  much  less  than 
the  plate) , 


*Some  additional  details  may  be  found  in  Reference  7. 
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It  is  evident  from  Fig.  1  that  doubling  of  the  distance 
between  screws  shifts  the  first  absorption  peak  by  a  factor 
of  4  in  frequency,  as  one  would  predict  on  the  basis  of  the 
aforementioned  resonances.  The  measured  values  of  the  first 
few  peak  frequencies  also  agree  surprisingly  well  with  values 
computed  as  if  the  plate  portions  between  bolts  were  simply 
supported  at  the  bolts.  As  has  been  remarked,  this  type  of 
calculation  a.ssumes  that  the  beam  acta  essentially  like  an 
infinitely  rigid  and  massive  support  for  the  plate  segments. 
The  good  agreement  observed  may  be  ascribed  to  the  fact  that 
the  beams  used  are  very  stiff  and  heavy  as  compared  to  the 
plate;  worse  agreement  should  be  expected  for  lighter  or 
softer  beams,  which  may  make  greater  contributions  to  tht* 
relative  motion. 

Parallel  beams,  with  dissipative  material  between  them  -  For 
systems  of  this  type,  shown  schematically  in  Figs.  1  £ind  2,  peak 
absorption  may  be  expected  to  occur  at  a  resonance  frequency  deter¬ 
mined  by  the  masses  of  the  beams  and  the  spring  action  of  the  plate 
between  them.  The  dissipative  mechanism  here  is  provided  by  defor¬ 
mation  of  the  damping  compound  (Aquaplas).  Frequencies  of  maximum 
absorption  calculated  on  this  basis  agree  well  with  the  observed 
values. 


One  may  observe  that  at  high  frequencies  the  glued-on 
beams  resulted  in  low  absorption  (Plg.  l),  whereas  the  bolted- 
on  beams  produced  fairly  high  absorption  (Fig.  2).  This  may 
be  due  to  the  glued-on  beams'  reflecting  most  of  the  incident 
bending  wave  energy,  which  thus  can  not  be  dissipated  by  the 
damping  material  between  the  beams.  On  the  other  hand,  flex¬ 
ural  waves  which  are  much  shorter  than  the  distance  between 
bolts  pass  the  bolted-on  beam  virtually  unattenuated  (ref.  8) 
and  can  therefore  dissipate  much  of  their  energy  in  the  damping 
material. 

Riveted  Beams  -  Measurements  indicated  in  Pig.  2  were  performed 
on  a  panel  made  to  approximate  aircraft  construction.  This  panel 
is  shown  in  Fig.  3.  Absorption  here  was  probably  due  to  friction 
betvfeen  the  beams  and  the  plates.  The  low  absorption  coefficients 
observed  might  be  due  to  the  close  rivet  spacing,  which  resti'icts 
the  relative  motions  of  beams  and  plate.  See  also  Section  Il-F.  of 
this  report. 
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B.  Structures  Excited  By  Reverberant  Acoustic  Fields; 

Application  of  Energy  Concepts _ _ _ 

Thermodynamic  or  energy  concepts  are  generally  useful  for 
arriving  at  some  tmderstanding  of  complex  systems,  since  these 
concepts  do  not  rely  on  detailed  analyses  of  the  systems  Involved. 
These  concepts,  properly  interpreted,  might  therefore  be  used  to 
advantage  in  analyzing  the  vibrations  of  complex  structures  ex¬ 
posed  to  random  acoustic  fields. 

The  average  energy  associated  with  a  given  mode  of  a  struc¬ 
ture  or  of  a  reverberant  room  may.  In  analogy  to  gas  dynamics,  be 
taken  as  a  measure  of  the  "acoustic  temperature"  of  the  structure 
or  room  in  a  small  frequency  Interval  enclosing  the  modal  frequency. 
One  may  then  conceive  of  a  reverberant  room  as  a  thermal  bath  in 
which  the  structure  is  immersed  and  with  which  it  will  eventually 
come  to  thermal  equilibrium.  The  equilibrium  temperature  of  the 
structure  then  depends  on  the  temperature  of  the  room  and  on  how 
good  the  thermal  contact  between  the  two  is.  This  thermal  analogy 
has  been  useful  in  guiding  the  analyses  presented  subsequently, 
but  further  elaboration  on  it  will  be  omitted  here  since  it  is  not 
essential  for  summary  and  application  of  the  results  of  these 
studies  and  since  it  is  available  in  detail  In  reference  9. 


1.  Basic  Concepts 

a)  Modal  Behavior 

vnienever  structural  resonances  occur  in  separate  sharply 
defined  frequency  regions;  i.e.,  whenever  the  effective  struc¬ 
tural  damping  is  very  small  (or  the  quality  factor  Q,  is  very 
high),  the  various  structural  modes  may  be  considered  as  inco¬ 
herent  harmonic  oscillators.  That  is,  each  mode  has  associated 
with  it  its  own  modal  velocity  energy,  mass,  force,  etc. 


In  complete  analogy  to  relations  developed  for  a  slngle- 
degree-of -freedom  system  subject  to  random  excitation,  one  may 
express  the  mean  aqxmre  modal  velocity  ^  as 


_  TT 

"2  M 


(B.l) 


where  3^(a))  denotes  the  power  spectral  density  of  the  modal  force 
f,  03  the  natural  frequency  associated  with  the  mode,  Rm  the  total 
resistance  (analogous  to  the  usual  viscous  damping  coefficient), 
and  M  the  modal  mass  (refs.  10-12), 
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For  structures  Immersed  In  acoustic  media  it  is  convenient 
to  consider  the  total  resistance  as  the  sum  of  two  terms t  one 
(R  ,)  to  account  for  energy  radiated  to  the  meditam,  and  the 
otnlr  to  account  for  other  energy  losses,  such  as  those 

due  to  material  hysteresis  or  "mechanical"  damping. 


Prom  the  work  of  Smith  (ref.  10 )  one  may  show  that  for  a 
structiire  exposed  to  a  reverberant  acoustic  field  with  power 
spectral  density  of  pressure  S  (co),  the  previously  defined  S„ 
is  given  by  ^ 


SfK) 


Uv  c 


2  ^rad 


P«  <*>1^ 
0 


(B.2) 


Ifere  p  denotes  the  density  of  the  acoustic  medlian,  c  the  speed 
of  sound  in  it. 

By  combining  Eqs.  (B.l)  and  (B.2)  one  may  obtain  the  follow¬ 
ing  cogent  expression; 

flg  =•  ®R  *  M-  (B-3) 

where 

^2 

Sq  =  M  V  =  average  energy  per  structural 

^  mode  (B.4) 

2Tr^  c 

0p  - - 5=  B  (o)  ni  average  energy  per  acoustic 

^  03^  P  mode*  (B.5) 

CL 


♦The  total  acoustic  energy  in  a  reverberant  room  of  volume  V  in 
a  frequency  Interval  Ao,  within  which  the  mean  square  pressure  is 

p  «Sp((o^)«Aa),  la  given  by 

Er  -  ^/pA 

and  the  number  of  modes  in  the  room  in  this  frequency  Interval 
is  (ref.  5) 

%  “o^  Aa3/2ir^  c|  . 

Bence  =»  given  in  Eq.  (B.5). 
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[j,  =  - Tu -  =  a  factor  that  denotes  the  degree  (B.6) 

^rad"*"  mech  of  coupling  between  room  and 
structural  mode. 

Thusj  If  the  various  tei'ins  that  enter  Eq.  {B.3)  are  In  hand  one 
may  utilize  this  equation  to  estimate  the  response  of  structures  to 
reverberant  excitation  under  the  conditions  where  the  assumptions 
underlying  the  derivation  of  Eq.  (B.3)  are  reasonably  well  satisfied. 
The  remainder  of  this  section  will  therefore  be  concerned  primarily 
with  evaluation  of  the  various  parameters. 

b)  Parameter  Estimates 
l)  Modal  Mass 

The  modal  mass  M  Is  defined  so  that  M  v  /2  gives  the  average 
kinetic  energy  of  the  structure.  The  modal  mass  Is  between  1/2  and 
1  times  the  total  mass,  it  is  nearer  l/2  of  the  total  mass  for  higher 
modes . 


ll)  Coupliiig  Factor  p. 

It  Is  evident  from  Eq.  (B.6)  that  for  those  cases  where  mechan¬ 
ical  damping  Rmech  is  much  smaller  than  radiation  damping  R-q^, 

For  most  practical  cases,  however,  ^rad'^'^\*ech^ 


_  ^rad 


(B.7) 


Equation  (B.3)  applies  when  an  acoustic  field  excites  the  structure. 
For  the  case  where  the  structure  excites  the  acoustic  field  one  finds 
from  a  relatively  detailed  examination  of  the  intei^actlon  at  e„=d_ix’ 
with  |i'  au.  ”  ^ 


ill)  Mechanical  Resistance 

As  for  asimple  harmonic  oscillator,  the  average  power  loss  for 
a  mode  is  P=Rq,^  and  the  average  energy  Is  E=Mv^.  By  equating  P/E 
obtained  here  to  the  same  ratio  computed  from  Eq.  (A. 4)  and  noting 
that  generally  Rinech»Bi,ad  express  the  mechanical  resistance 

directly  In  terms  of  the  previously  introduced  structiiral  loss  factor 
T|  and  absorption  coefficient  7: 

Voh"  ^  1' ]  •  tB.8) 
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ir  boundary  losses  predominate  over  internal  losses,  as  often  is  the 
case,  then 


R  «  ffi 

Wch  TT 


M 

S 


(B.9) 


where  7^  is  the  absorption  coefficient  of  a  portion  of  the  boundary 
of  length  LjL  and  M/S«ps/2  (i.e.,  M/S  is  approximately  equal  to  one 
half  the  structural  mass  per  unit  surface  area) . 


iv)  Pladiation  Resistance 

The  coincidence  phenomenon  is  known  to  be  a  major  factor  in 
determining  the  radiation  resistance  of  simple  structures.  At 
frequencies  below  the  coincidence  frequency  the  length  of  a  flex¬ 
ural  wave  in  a  structure  is  greater  than  the  acoustic  wavelength  in 
the  surrounding  medium,  and  the  coupling  between  the  structure  and 
the  medium  is  poor;  RT>ad  small.  At  frequencies  above  coincidence 
the  flexural  wavelengths  are  smaller  than  the  acoustic  wavelengths, 
coupling  is  better,  and  is  greater. 

For  imiform  beams  or  plates  the  coincidence  frequency  f^  is 
given  by 

2irf^  =  “c  ~  '  (B.IO) 

where  k  denotes  the  radius  of  gyration  of  the  cross-section 
(k  =  h/ij  12  for  plates  of  thickness  h) ,  c  denotes  the  speed  of 
sound  waves  in  the  acoustic  medium,  and  cr  the  speed  of  longitud¬ 
inal  waves  in  the  structural  material.  Stiff er  structural  members 
thus  have  lower  coincidence  frequencies  and  might  be  expected  to 
exhibit  higher  average  radiation  resistances  in  a  given  wide  fre¬ 
quency  band.  •* 

In  a  structure  composed  of  many  beams  and  panels  (e.g.,  a 
typical  aircraft  section)  one  may  generally  expect  to  find  that 
some  components  of  the  structure  are  below  and  others  above  coin¬ 
cidence  at  a  given  frequency.  In  order  to  obtain  a  first  approx¬ 
imation  to  the  radiation  resistance  of  such  a  composite  structure 
one  might  neglect  interaction  effects  and  assume  that  only  those 
components  which  are  above  coincidence  contribute  significantly  to 
the  radiation  resistance.  Radiation  resistance  estimation  proce¬ 
dures  based  on  these  assumptions  are  outlined  below. 


Since  the  beams  in  most  practical  beam-plate  structures  are 
much  stiffer  than  the  plates,  the  beams  should  dominate  the  radia¬ 
tion  resistance  of  such  structures  at  the  lower  frequencies.  The 
panels,  by  virtue  of  their  higher  coincidence  frequencies,  may  have 
little  effect  except  at  the  highest  frequencies. 
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The  aco’ustlc  radiation  of  beam-plate  systems  in  which  the  plates 
are  limp  (i.e.,  soft  and  damped)  may  thus  be  assumed  to  be  primarily 
due  to  the  beams  acting  as  strip  radiators,  with  the  panels  acting 
essentially  as  baffles,  but  perhaps  adding  a  little  to  the  width  of 
the  strip  radiators  -  in  analogy  to  a  similar  situation  encountered 
in  studies  of  sound  transmission  through  walls  with  studs  (ref.  13)* 

The  radiation  resistance  of  a  strip  radiator  of  length  L  and 
width  w  is  (for  both  faces)  given  by  (ref,  l4) 

"strip  “  '  Poi-  “  >  “o 

provided  that  the  strip  width  w  is  small  compared  to  the  acoustic 
wavelength.  Below  coincidence  (that  is,  for  c><(s  )  the  radiation 
resistance  of  a  strip  is  much  smaller  than  that  given  above. 

Heckl  (ref.  15)  has  suggested  that  a  strip  radiator  attached  to 
a  limp  plates  has  an  effective  radiating  width  wwXp/ir,  where  A-  is 
the  wavelength  of  plate  flexural  waves.  For  homogeneous  plates, 

Xp  =  Sir  Kc^oi  (B.I2) 

so  that  for  a  beam  attached  to  limp  panels, 

^strip  “  (B.13) 

where  k  and  Ct  refer  to  plate  radius  of  gyration  and  longitudinal 
wave  velocity. 

One  may  note  that  Tvp  decreases  with  increasing  frequency  and 
that  Hg-tpjLp  is  Independent  of  frequency.  However,  the  effective 
strip  width  associated  with  a  practical  beam  probably  will  not  be 
less  than  the  width  of  the  beam  flange  in  contact  with  the  plate. 

If  this  flange  width  is  denoted  by  w  ,  then  ^  frequency 

fp  given  by  o  p' 


^r  = 


2kc, 


TM, 


(b.i4) 


For  f>fp  the  radiation  resistance  of  the  strip  should  obey 

"atrip  “  “LPa-”o  ' 
instead  of  Eq.  (B.13)» 


(B.15) 
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For  a  beam-plate  structure  In  which  the  plates  are  below 
coincidence  the  radiation  resistance  at  frequency  f  may  then  be 
approximated  by 


where  the  first  summation  denotes  the  total  length  of  beams  that  at 
frequency  f  are  above  coincidence  but  below  f  ,  and  the  second  sum¬ 
mation  la  taken  over  all  beams  for  which  f>f^. 

On  the  other  handj  at  frequencies  which  are  above  th^  coincidence 
frequency  of  the  panels,  one  may  expect  the  radiation  to  be  dominated 
by  that  due  to  the  panel.  Then,  If  one  neglects  edge  effects  (ref.  l) 
and  takes  Into  account  radiation  from  both  faces,  where  A  Is  the  total 
area  of  the  panels  (one  face)  that  are  above  coincidence, 

•Vad  ”  ,  (B.17) 


It  must  be  pointed  out  that  In  computing  radiateiLpower  from  the 
defining  relation  ^  one  must  always  use  that  v^  which  refers 

to  the  velocity  of  the  radiating  section.  Since  It  la  the  panel  that 
radiates.  Its  velocity  must  be  used  In  the  foregoing  expression. 

(When  the  panel  Is  "limp"  the  velocity  of  that  portion  of  It  that  Is 
radiation  Is  essentially  that  of  the  beam.  Validity  of  the  panel 
velocity  under  this  condition  has  been  demonstrated  eiqperlmentally, 
as  discussed  In  the  next  section.) 


A  relation  between  panel  velocity  Vp  and  beam  velocity  v^  may  be 

deduced  for  the  condition  of  modal  Independence  from  the  assui^tion 
that  the  average  energy  per  mode  is  the  same  for  all  stmictural  com¬ 
ponents.  Then  the  energy  contributed  by  a  component  is  proportional 
to  its  modal  density,  and  If  the  beams  and  the  plate  are  of  the  same 
material. 


“bi 


(B.18) 


where  Mp  denotes  the  total  panel  mass,  the  total  beam  mass,  and 

denotes  the  plate  radius  of  gyration,  the  beam  radius  of  gyration. 
(This  equation  assumes  the  grillage  composed  of  the  beams  to  behave 
like  an  equivalent  plate  (ref.  16; .  The  atlffest  beams  will  Influence 

this  behavior  most,  therefore  their  k  values  should  be  used  In  this 

P 

approximate  computation.)  If  beams  and  plates  are  of  different  mate¬ 
rials,  the  right-hand  side  of  Eq.  (B.18)  should  be  multiplied  by 
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It  should  be  remembered  that  the  foregoing  procedure  for  estimat¬ 
ing  radiation  resistance  constitutes  only  of  a  first  step,  and  as  such 
may  be  unreliable.  Pertinent  results  of  a  later  study  (ref.  17)  are 
discussed  under  the  heading  of  Conclusions  at  the  end  of  this  section 
(page  22). 

c)  Non-Modal  Behavior 

i)  Limit  of  Modal  Behavior 

The  entire  previous  discussion  is  based  on  the  assumption  that 
the  structural  modes  do  not  overlap  In  frequency,  so  that  their  vibra¬ 
tions  are  statistically  independent  and  their  energies  do  not  Interact. 
Overlap  occurs  If  the  frequency  bandwidth  associated  with  a  mode  Is  of 
the  same  order  of  magnitude  as  the  average  frequency  separation  between 
modes. 

As  shown  in  Appendix  I,  the  plate-like  portions  of  beam-plate  sys¬ 
tems  contribute  many  more  modes  than  the  beam-llke  portions,  so  that 
the  former  predominate  in  establishing  the  separation  between  modes. 

The  average  frequency  separation  Af  between  modes  of  a  panel  of  area  A 
and  thickness  h  is  given  by 

Af  «  hCj/Ayi”*  ,  (B.I9) 

where  Cj^  denotes  the  longitudinal  wave  velocity  In  the  panel  material. 
The  same  expression  applies  also  for  the  higher  frequencies  of  a  beam- 
plate  system  in  which  all  panels  are  of  the  same  thickness  and  material; 
A  then  denotes  the  total  panel  area. 

If  one  assumes  for  the  sake  of  simplicity  that  the  location  of 
structural  modes  along  a  frequency  axis  may  be  described  by  a  Poisson 
process  (ref.  I8),  one  may  show  that  the  probability  of  finding  more 
thaux  one  mode  in  a  frequency  interval  Ab  Is  0.5  or  greater  if 
Ab/Af>1.68.  One  may  then  define  a  "frequency  of  modal  overleip"  fco 
a  structure  as  that  frequency  for  which  the  modal  bandwidth  equals 
1.68  Af.  Since  bandwidth  is  related  to  quality  factor  ft  or  loss  factor 
tj  as  Ab»f/ft=*xif,  the  frequency  of  modal  overlap  is  given  by 

f^  «  1.68ftAf  =  1.68  ^  (B.20) 

At  frequencies  above  this  frequency  the  probability  of  finding  more 
than  one  modal  peak  within  bandwidth  Ab  of  a  modal  response  is  greater 
than  0.5.  Thus,  modal  behavior  of  the  structure  may  reasonably  be 
assumed  at  frequencies  below  f^. 
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11 )  Estimation  of  Non-Modal  Response 

At  frequencies  above  f  a  bejim-plate  structure  may  be  expected  to 
respond  to  an  acoustic  field  essentially  like  an  Infinite  plate,  and 
to  obey  the  "mass  law"  relation 

S^(cj)  =  4Sp(a))/p|  (B.21) 

where  S  is  the  power  spectral  density  of  acceleration  of  the  structure 
Pg  Is  Its  mean  mass  per  unit  surface  area. 

2.  Experimental  Results 

A  series  of  experiments  was  tindertaken  In  order  to  obtain  an  in¬ 
dication  of  the  validity  of  the  previously  outlined  theories.  These 
ejqjeriments  were  performed  primarily  on  a  simulated  aircraft  panel 
(shown  In  Pig.  3)  suspended  in  a  large  test  room.  In  each  case  the  ex¬ 
periments  were  performed  with  the  panel  In  Its  original  condition  and 
also  with  some  damping  (several  randomly  oriented  strips  of  MMH  damping 
tape)  added.  All  experiments  were  performed  with  narrow  band  noise  ex¬ 
citation. 

Figure  4  stimmarizes  results  obtained  from  experiments  consisting 
of  driving  the  panel  at  various  points  by  means  of  a  mechanical  shaker 
and  measuring  the  power  F  .  radiated  to  the  room  as  well  as  the  ve¬ 
locity  of  one  point  on  the“panel.  (A  point  on  a  rib  was  selected  since 
It  was  found  that  the  transducer  distorted  the  behavior  of  the  panel  at 
higher  frequency  when  attached  to  the  panel  itself.)  The  measured  quan 
titles  are  interpreted  In  the  figure  in  tarns  of  radiation  resistance 

^rad  relation  ^pad^^rad'^'^b*  subscript  b  here  in¬ 

dicates  that  velocity  was  measured  on  a  beam,  rather  than  on  a  panel, 
as  previously  indicated. 

The  results  of  reverberation  time  measurements  on  the  panel  of 
Fig.  3  are  summarized  In  Fig.  5,  interoreted  In  terms  of  total  resist¬ 
ance  ^*^inech'*"^rad  13. 81^^  with  the  modal  mass  M  taken  as 

1/2  the  total  mass. 

The  upper  series  of  points  shown  In  Pig,  6  represent  values  of  the 
coupling  factor  M-'  computed  according  to  Eq.  (B.6)  from  the  data  of 
Pigs.  4  and  5.  These  points  pertain  to  the  case  where  the  panel  is 
driven  mechanically  and  excites  acoustic  oscillations  In  the  test  room. 
The  lower  series  of  points  shown  In  Fig.  6  represents  values  of  the 
coupling  factor  p.  obtained  fTOm  additional  measurements  In  which  sound 
In  the  room  excited  vibrations  In  the  panel.  The  latter  results  were 
computed  on  the  basis  of  modal  energies,  using  Eqs.  {B.3)  and  (B.5). 
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A  DRIVEN  AT  POINT  8  OF  FlG.  3 

A  DRIVEN  AT  POINT  4  OF  FIG.  3 

O  DRIVEN  AT  POINT  4  OF  FIG.  3 
WITH  ADDED  DAMPING 
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FIG.  4  RADIATION  RESISTANCE  OF  SIMULATED  AIRCRAFT  PANEL,  AS 
DETERMINED  FROM  MEASUREMENTS  OF  POWER  RADIATED  TO 
ROOM  AND  VELOCITY  OF  POINT  6  OF  FIG.  3 
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COUPLING  FACTORS  ft  AND^'  FOR  SIMULATED  AIRCRAFT 
PANEL  (Vjj  MEASURED  AT  PT.6  USED  FOR  PANEL  VELOCITY) 


All  computations  for  Fig.  6  were  based  on  the  assumption  that  as 
measured  at  point  6  shovm  In  Fig.  3,  represents  the  velocity  of “the 
simulated  aircraft  panel.  Theory  Indicates  that  p.  and  p'  of  Fig.  6 
should  be  equal 3  and  it  appears  that  panel  velocity  Vp  Instead  of 
beam  velocity  v^  should  be  used  in  the  various  computations.  The  ra¬ 
diation  resistance  referred  to  a  beam  is  related  to  that  referred 

to  a  plate  according  to 


''p 


(B.22) 


Therefore  the  p*  data  of  Fig.  6  appears  to  be  too  high  by  a  factor  of 

v^/v^  and  the  p  data  appear  to  be  too  low  oy  the  same  factor.  Hence 

this  factor  may  be  computed  by  comparison  of  the  two  sets  of  data.  The 
results  of  such  a  computation  are  shown  in  Fig.  Tj  together  with  the 
value  estimated  from  Eq.  (B.I8)  which  is  based  on  the  assumption  of 
equipartltion  of  energy  among  all  structural  modes. 

Figure  8  shows  the  data  of  Pig.  4  corrected  from  to  by 

use  of  the  velocity  ratio  plotted  In  Pig.  7.  Also  shown  in  Pig,  8  are 
the  radiation  resistance  values  predicted  by  Eq.  (B.16). 


3.  Conclusions 


Prom  the  foregoing  experimental  results,  and  from  addition  results 
obtained  under  NASA  sponsorship  (ref.  I7)  subsequent  to  completion  of 
the  work  described  here,  one  may  arrive  at  the  following  conclusions: 

1.  Equipartltion  of  energy  among  the  structural  modes  is  reason¬ 
able  and  provides  a  useful  estimate  of  the  beam/plate  ve¬ 
locity  ratio. 

2.  The  modal  energy  expressions  of  Eqs.  (B.4)  and  (B.5)  appear 
to  be  valid,  and  the  modal  mass  may  reasonably  be  approxi¬ 
mated  by  1/2  the  total  mass. 

3.  Velocity  and  radiation  resistance  should  be  referred  to  the 
plate,  not  to  ribs,  in  calculation  of  coupling  factor  p.. 

4.  Reciprocity  holds;  l.e.,  p.’  (room  driving  plate)  =  n  (plate 
driving  room) . 

5.  Radiation  resistance  values  estimated  from  Eq.  (B.I6)  always 
exceed  the  actual  values,  occasionally  by  considerable 
amoxmts . 
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FIG.  7  RATIO  OF  MEAN  SQUARE  VELOCITIES  OF  PANELS  AND  BEAMS 
REQUIRED  TO  MAKE  DATA  OF  FIG.  6  SATISFY 
RECIPROCITY  (ii  =  u') 
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FIG.  8  RADIATION  RESISTANCE  VALUES  OF  FIG.  4,  CORRECTED  FOR 
PANEL/BEAM  VELOCITY  RATIO  OF  FIG.  7 


The  deviations  of  experimentally  detemlned  from  theoretically 
predicted  results  apparent  in  Fig.  8,  and  more  pronounced  deviations 
observed  In  a  later  study  (ref.  1?)  performed  with  a  larger  and  thicker 
baffled  plate  vrlth  more  massive  beams ;  have  shed  some  doubt  In  the  gen¬ 
eral  applicability  of  the  strip  radiator  concept  as  outlined  here.  A 
later  (so  far  only  partially  validated)  theory  (ref.  17)  ascribes  the 
radiation  resistance  of  beam-plate  structures  to  the  scattering  of 
plate  flexural  waves  at  beams,  edges,  or  other  discontinuities.  The 
presence  of  beams,  etc.,  thus  is  still  held  responsible  for  increas¬ 
ing  the  radiation  resistance  of  the  composite  structure,  but  beam  co¬ 
incidences  now  assume  much  less  importance.  Tlae  newer  model  predicts 
radii ation  resistance  increases  due  to  besun  coincidences  to  be  notice¬ 
able  only  v.’ith  beams  that  are  very  light  relative  to  the  plate  and  to 
occur  at  frequencies  which  may  be  considerably  higher  than  the  classi¬ 
cal  coincidence  frequencies  of  the  beams.  However,  since  the  strip 
radiator  theory  proposed  in  the  present  report  may  be  shown  on  the 
basis  of  the  newer  theory  always  to  result  in  high  estimates  of  the 
radiation  resistance,  the  simpler  strip  radiator  theory  may  still  be 
useful  for  conservative  design  estimates. 
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SECTION  II 


SPECIFIC  VIBRATION  CONTROL  TECHNIQUES 


The  maxlmvun  responses  of  structui'^es  to  vibratory  excitation  are 
affected  primarily  by  the  magnitude  of  the  excitation,  its  coupling  to 
the  structure,  and  by  the  resonant  and  damping  characteristics  of  the 
structure.  Control  of  vibration  must  therefore  be  concerned  with  fa¬ 
vorable  ad;}ustraent  of  these  factors . 

The  structural  designer  cannot  usually  hope  to  do  much  about  con¬ 
trolling  the  magnitudes  of  excitation  caused  by  aerodynamic  or  acoustic 
phenomena.  However,  when  the  excitation  is  a  local  one,  as  for  example 
due  to  a  vibrating  machine  (that  cannot  be  modified  itself),  the  de¬ 
signer  may  turn  to  vibration  absorbers  to  reduce  the  excitation  that 
reaches  the  structure.  Attempts  at  devising  vibration  absorbers  with 
improved  characteristics  are  summarized  in  the  first  of  the  following 
sections. 

Means  for  reducing  the  coupling  between  sources  of  excitation  and 
structures  are  discussed  in  the  next  two  sections.  Section  B  outlines 
experimental  results  that  demonstrate  the  utility  of  "velocity- controlled” 
beams  (whose  stiffness  decreases  with  increasing  frequency)  in  reducing 
the  radiation  resistance  of  beam-plate  structures.  Section  c  summa¬ 
rizes  results  of  a  mathematical  study  of  viscoelastic  leaf  spring  isola¬ 
tors,  vihlch  may  possess  some  attractive  features  for  reducing  the  cou¬ 
pling  between  local  sources  and  structures. 

The  final  three  sections  deal  with  aspects  of  structural  danqplng, 
which  is  most  Inqjortant  in  limiting  structural  responses.  Section  D 
outlines  a  theoretical  study  of  the  Interaction  of  beams  and  plates  in 
beam-plate  systems,  in  order  to  determine  how  application  of  dancing 
to  one  or  the  other  of  these  components  affects  the  system  responses. 
Section  E  presents  the  results  of  a  primarily  experimental  Investigation 
Intended  as  a  first  step  toward  understanding  and  improving  the  daaqjlng 
of  aircraft  structural  Joints.  Finally,  section  P  summarizes  the  deriva¬ 
tions  of  some  expressions  that  relate  the  loss- factors  of  beam-like 
structures  to  absorption  coefficients  exhibited  by  these  structures 
attached  to  plates. 

A.  Non-classical  Vibration  Abnorbers 

It  is  well  known  that  one  may  reduce  the  steady-state  response  of 
a  primary  mass  M  to  sinusoidal  excitation  at  frequency  co.  by  connecting 

^  2 

M  to  a  secondary  mass  m  via  a  spring  of  stiffness  k,  so  that  =  k/m. 

However,  the  secondary  mass  and  spring  can  reduce  (i.e.,  "absorb")  the 
vibration  of  the  primary  mass  only  for  frequencies  near  <0^.  For  exaitple. 
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for  driving  frequencies  only  slightly  greater  than  the  "absorber" 
results  in  large  amplification  of  the  response  of  M. 

1 )  Frequcincy-Variable  Springs 

If  the  frequency  of  absorption  could  be  made  to  change 

with  driving  frequency  co  so  as  to  remain  equal  or  nearly  equal  to  Itj 
the  resulting  device  would  act  Hkp  an  ideal  absorber  over  the  entire 
range  over  which  this  equality  is  maintained. 

Continuous  reduction  of  the  secondary  mass  with  increasing  fre¬ 
quency  appears  impossible  to  obrain  practically.  However,  viscoelastic 
springs  are  known  to  increase  in  stiffness  with  increasing  frequency. 
Detailed  consideration  of  the  absorber  action  of  viscoelastic  springs 
6ind  of  representative  suitable  materials  showed  that; 

1)  The  stiffnesses  of  viscoelastic  materials  Increase  at 
most  as  the  first  power  of  frequency,  whereas  Increase 
as  the  square  is  required  for  continuous  optimum  ab¬ 
sorption. 

2)  The  most  rapid  increases  of  viscoelastic  material  stiff¬ 
ness  occur  in  conjunction  with  high  damping,  which  de¬ 
creases  absorber  effectiveness. 

3)  Large  ratios  of  secondary  to  primary  mass  would  be  re¬ 
quired  for  response  reduction  over  a  wide  frequency 
range.  Large  added  masses  are  generally  undesirable; 
otherwise  they  could  be  connected  rigidly  to  the  primary 
mass  to  reduce  Its  response  without  introducing  new 
resonances. 

These  conclusions,  coupled  with  practical  considerations  such  as  the 
temperature  dependence  of  the  properties  of  viscoelastic  materials, 
tend  to  rule  out  the  desirability  of  vibration  absorbers  with  visco¬ 
elastic  springs. 

2)  Distributed  Secondary  System 

It  was  thought  that  an  absorber  using  a  distributed  mechanical 
system  (such  as  a  beam)  instead  of  a  rigid  secondary  mass  m  could  be 
useful  for  controlling  the  vibrations  of  the  primary  mass  at  a  nxanber 
of  frequencies.  Detailed  analyses  of  this  idea  showed  that  absorption 
could  Indeed  be  attained  at  many  frequencies,  but  that  such  absorbers 
also  result  in  resonances  (and  severe  amplification  of  the  responses) 
at  frequencies  very  near  the  absorption  frequencies.  Such  systems 
would  therefore  require  very  precise  absorber  design  and  accurate  pre¬ 
diction  of  the  exciting  frequencies,  and  thus  must  be  Judged  to  be  im¬ 
practical  in  general. 
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B •  Had latlon  Resistance  Control  in  Be am  Design* 

Beams  may  be  designed  so  that  their  flexural  rigidities  de¬ 
crease  with  increasing  frequency.  They  thus  have  higher  coincidence 
frequencies  than  conventional  beams  with  the  same  low-frequency  (static) 
stiffness  and  should  therefore  result  in  less  coupling  to  an  acoustic 
field. 

Two  4^1"  X  24"  test  panels  of  0.032"  thick  aluminum  v;ere  construct¬ 
ed  for  this  purpose.  Three  24"  long  aluminum  beams  with  C  shaped 
cross-sections  (O.O65"  thick  1,',''5”  deep  1.0"  flanges)  were  attached 
to  one  of  the  panels  by  means  of  an  epoxy  adhesive  parallel  to  the 
24"  panel  dlinetialon  and  uniformly  spaced  over  the  44"  length.  Three 
beams  designed  as  indicated  in  Reference  19  of  approximately  the  same 
static  stiffness  as  the  aluminum  channels  were  similarly  attached  to 
the  second  panel.  The  coincidence  frequency  of  the  aluminum  channels 
was  calculated  to  be  about  200  cps  that  of  the  latter  type  beams 
about  5  ’■<c . 

The  results  of  measurements  carried  out  on  these  two  test  panels 
appear  in  Fig.  9.  The  latter  beams  were  found  to  produce  lower 
values  of  radiation  resistance  and  coupling  factor^  than  the 

more  conventional  channels.  The  near  equality  of  the  total  resist¬ 
ances  H^q^  of  the  two  panels  indicates  that  the  mechanical  resistances 

R  h  were  nearly  equal  and  considerably  greater  than  the 

respective  radiation  resistances. 

C ,  Viscoelastic  Vibration  Isolators 

1.  Leaf  Springs 

Viscoelastically  damped  beams  may  be  designed  according  to 
available  methods  (refs.  2  20)  to  satisfy  a  considerable  range  of 
requirements  and  thus  might  provide  superior  broadband  vibration 
isolation  characteristics.  An  exploratory  study  of  viscoelastic  leaf 
springs  vjas  undertaken,  therefore,  as  discussed  in  the  following  pages. 

Consider  a  mass  mounted  symmetrically  on  tv;o  identical  leaf 
springs.  If  a  harmonic  force  of  amplitude  F  acts  on  this  mass,  then 
a  (rigid)  support  to  which  the  springs  are  attached  will  also  experience 
a  harmonic  force,  but  of  amplitude  F^.  Alternatively  if  the  support 
is  made  to  oscillate  harmonically  vjirh  an  displacement  amplitude  Yq. 

then  the  supported  mass  will  oscillate  with  an  amplitude  Y.  The  trans- 
missibility  of  a  mounting  system  (here  the  leaf  springs)  is  defined  as 


*The  experiments  described  here  use  beams  which  are  the  subject  of  pat¬ 
ent  application  by  Bolt  Beranek  and  Newman,  Inc,,  and  which  bear 
the  trade-mark  "Soundshear" .  Properties  and  design  of  these 
structures  are  described  in  Reference  19. 
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the  equivalence  of  the  two  ratios  may  be  deduced  from  reciprocity 
(ref.  21). 

The  characteristics  of  beam-like  structures  depend  on  their  bound¬ 
ary  conditions.  For  the  present  analyses  the  leaf  springs  were  as¬ 
sumed  to  be  ideally  clamped  at  the  support  and  at  the  mounted  mass. 

By  solving  the  equation  of  flexural  motions  of  elastic  beams  (refs.  4, 
21)  subject  to  the  appropriate  boundary  conditions  one  finds  that 
(ref,  22) 


T  = 


h(a) 

f(a)  +  li*g(ct) 


where 


(C.2) 


h(a)  *  sinh  a  +  sin  a 

f(a)  “  cos  a* sinh  a  +  sin  u-cosii  a 

(C.3) 

g(a)  «  a(coB  a  cosh  a  -  l) 
a  *  kL  =  2ir-L/^ 

M.  =  M/mL 

Here  L  denotes  the  length  of  one  beam,  m  its  mass  per  unit  length,  and 
M  the  supported  mass  per  beam  (or  half  the  supported  mass).  Thus,  u, 
denotes  the  ratio  of  supported  mass  to  the  mass  of  the  supporting  de¬ 
vice..  The  symbol  k  denotes  the  wave  number,  7^  the  wavelength  associa¬ 
ted  with  free  flexural  waves  in  the  beam; 


k 


mtc^EI 


E  denotes  Young's  modulus  of  the  beam  material, 
of  Inertia  of  the  beam  cross-section. 


(C.4) 

I  the  oentroldal  moment 


One  may  take  damping  into  acco\ant  by  letting  E  take  on  con^jlex 
values;  that  is,  by  replacing  E  by  E*  =  +  lE^  =  E^  (1  +  1^)  where 

Eg  denotes  the  loss  modulus,  the  loss  factor  of  the  beam  (refs.  2, 

23).  k  and  a  also  take  on  complex  values  k*  and  a*. 

The  general  algebraic  expressions  tend  to  become  intractable,  but 
fortunately  the  mass  ratios  p.  encountered  in  practical  situations  tend 
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to  fall  into  two  categories;  l)  ix  «  1  (mass  of  leaf  spring  pre¬ 
dominates)  and  2)  u  »  1  (supported  mass  predominates).  For  these 
extreme  situations  one  may  arrive  at  some  relatively  simple  approx¬ 
imations. 


For  p,  »  1  and  small  a  (low  frequencies)  one  finds  that  Eqs. 
fC.2)  and  (C.3)  reduce  to  the  classical  transmissiblllty  expression 
(refs,  22,  24)  whose  behavior  is  sketched  in  the  Insert  of  Fig.  10. 
For  large  a  and  large  p,  one  obtains  the  appi’oxiraatlon 


a*  P.T 


sin  g*  +  sinh  a* 
cos  a*  cosh  a*  -  1 


A _ 

j^s  a* 


(C.5) 


The  second  appi’oximate  equality  applies  for  |a*  |  not  ^oo  small,  -  a 
condition  which  may  be  verified  to  apply  where  wave  effects  are  im¬ 
portant.  If  one  introduces  the  parameters  A  and  0,  defined  so  that 


A  =  |a*|  ,  0  =  ^  tan  q  ,  a*  ==  a'  -  la"  =  Ae~^^  (C.6) 


then  one  finds  that  transmissiblllty  peaks  (standing  wave  resonances) 
occur  where  Re£cos  =  cos  a'  cosh  a"  =  0,  or  for 

a'  =  7r(n  +  ^)  ,  n  =  1,  2,  ♦  •  •  .  (C.7) 

The  magnitudes  T  of  the  transmissiblllty  peaks  are  then  given  by 

p,^*  «  (A  sinh  a")"^  «  [a^  sinh  (a^  sin  0)]“^  .  (C.8) 

For  most  practical  calculations  one  may  determine  the  transmissiblllty 
peaks  from  Eqs.  (C,7)  and  (C.8),  In  addition  A»a'  unless  the  damping 
is  extremely  high*,  and  practical  calculations  may  be  simplified  fur¬ 
ther. 


The  validity  of  the  approximations  of  Eqs.  (C.7)  and  (C,8)  may 
be  verified  from  Fig.  10,  in  which  the  curves  represent  the  results 
of  calculations  carried  out  by  means  of  a  digital  computer  (ref.  25) 
on  the  basis  of  the  more  exact  expression  Eq.  (C.5)  and  the  points 
marked  by  crosses  represent  values  calculated  from  the  approximations 
Eqs.  (c.7)  and  (c,8) 


For  p,  «  1  and  a  not  too  small  one  similarly  may  write 


T  « 


_ ^ _ 

"  jsin  a*  +  cos  a* 


(C.9) 


This  function  exhibits  peaks  approximately  where  the  real  part  of  the 
denominator  vanishes,  or  for 


*For  q  ^  2,0  one  finds  1  ^  a '/A  1.04. 
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a*  =  Tr(n  + -J)  ,  n  =  0,  1,  2,  (C,10) 

A 

The  transmlssibllity  peak  values  T  may  be  found  by  substituting  the 
foregoing  values  of  a'  into 

T  «  (2  sin  ex'  sinh  a")”^  »  [ -^  slnh{a^  sin  <#>)]~^,  (C.ll) 

Figure  11  shows  curves  calculated  fi’om  th?-  more  exact  expression 
Eq.  (C.9)  and  peak  value  points  obtained  from  the  foregoing  approx¬ 
imations.  Good  agreement  is  again  evident. 

2.  Axial  Springs 

It  is  instructive  to  compare  the  results  obtained  for  leaf 
springs  with  similar  results  for  axial  (e.g.,  compression)  springs. 
For  the  latter  the  general  transmisslbility  expression  may  be  written 
as  (refs.  22,  26} 

T  =  I  cos  p  -  sin 


U'  =  M/mH 


where  m  denotes  the  mass  per  unit  length,  H  the  equilibrium  length 
of  the  spring,  K  the  classical  spring  constant,  and  n'  the  ratio  of 
supported  mass  to  total  spring  mass. 

One  may  again  introduce  dsimping  effects  by  replacing  K  by  an 
appropriate  complex  quantity  K*  =  K,  +  1K„  =  K^Cl+lti);  then  p  also 
takes  on  complex  values;  ^ 

p*  =  p'  -  ip"  =  Be"^'^  ,  B  =  1p*1  ,  7  =  "I  arc  tan  q 

(C.13) 

One  again  finds  that  for  large  mass  ratio  u'  and  low  frequen¬ 
cies  (small  p*)  Eqs.  (G.12)  reduce  to  the  classical  transmisslbility 
result  (ref.  22).  For  p*  not  too  small  and  large  p,  Eq.  (C.12)  may 
be  replaced  by  the  approximate  relation 
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mL=  BEAM  MASS,  M=MOUNTED  MASS,  K  "•  STATIC  STIFFNESS,  ii  =  K/M,  ;i=M/mL 

FIG.lt  TRANSMISSIBILITY  OF  VISCOELASTIC  LEAF  SPRING  FOR 
SMALL  MOUNTED  MASS  l.O) 


Cc.14) 


Ip'^Ih'T  s:  I  Sin  P*l"^ 

which  implies  transmiaslbllity  peaks  given  by 

H'T  =  (B  sinh  p")"^  «  [B  sinh  (B  sin  7)]“^  (C.I5) 

whene 


P'  s:  rnr 


n  =  1, 


00* 


Similarly;,  for  small  mass  ratio  p,', 


{ G , 16 ) 


T  ~  |cos  p*|"^  ,  (C.17) 

T  w  I sinh  p"l”^  =  [sinh  (B  sin  7)]“^ 
and  peaks  occur  where 

P '  «  Tr(  n  +  ^)  ,  n  =  1,  2,  •  •  •  { C .  I8 ) 

Tiie  validity  of  these  approximations  may  be  inferred  from  Pigs, 
12  and  13,  in  which  the  curves  represent  the  results  of  the  more 
exact  calculations  and  the  points  represent  results  obtained  by  means 
of  the  foregoing  approximate  expressions. 


3 .  Comparison 

It  may  be  shown  that  if  leaf  springs  and  compression  springs  are 
to  have  the  same  ratio  of  supported  mass  to  support  mass  and  the  same 
classical  fxnndamental  natiiral  frequency,  then  the  previously  defined 
parameters  A  and  B  must  be  related  as 


Figures  11  and  12  are  scaled  to  conform  to  Figs,  9  and  10  on  this 
basis,  so  that  one  may  compai'-e  leaf  springs  and  common  springs 
equitably  by  comparing  these  sets  of  figures. 


On  the  basis  of  these  figures,  or  by  means  of  mathematical 
analysis  of  the  pertinent  relations,  one  may  conclude  that  compared 
to  compression  springs  of  equivalent  mass  and  damping,  leaf  springs 

1)  result  in  fewer  transmisslblllty  peaks  for  a  given 
frequency  range ;  but 

2)  tend  to  result  In  a  generally  higher  level  of  trans- 
missiblllty  between  peaks  and  In  a  higher  envelope  of 
transmlsslbility  peaks. 
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(MASS  RATIO)  ■  (TRANSMISSIBILITY) 


FiG.l2  TRANSMISSIBILITY  OF  VISCOELASTIC  COMPRESSION 
SPRING  FOR  LARGE  RATIO  /i' OF  MOUNTED 
MASS  TO  SPRING  MASS 


6 
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ELASTIC 
MASS  C/i' 


One  should  note,  hiowever,  that  the  present  state  of  the  art 
permits  large  amounts  of  damping  to  be  designed  into  leaf  springs 
more  readily  than  into  compression  springs,  and  that  viscoelastic- 
ally  damped  leaf  springs  may  thus  constitute  superior  broadband 
vibration  isolation  devices. 

Damping  of  Beam-Plate  Systems 

Since  aircraft  and  other  structures  consist  primarily  of  beams 
and  plates,  the  reduction  of  the  vibrations  of  such  beam-plate  sys¬ 
tems  has  received  considerable  attention.  The  addition  of  damping 
structures  and  materials  usually  results  in  weight  Increases,  so  that 
it  is  desirable  to  apply  such  damping  treatments  as  conservatively  as 
possible  and  only  where  they  will  do  the  most  good.  The  subsequently 
outlined  study  was  undertaken,  therefore,  as  an  attempt  to  define 
under  what  conditions  damping  of  only  the  beams  is  preferable  to  damp¬ 
ing  of  only  the  plates,  and  vice  versa. 

The  system  studied  is  sketched  in  Fig,  l4.  Although  this  sketch 
shows  a  rectangular  plate  and  the  subsequent  discussion  deals  in 
Cartesian  coordinates,  these  restrictions  are  not  necessary;  the 
conclusions  of  the  study  should  be  more  generally  applicable.  However, 
it  must  be  assumed  that  the  boundary  conditions  of  the  beam  are  the 
same  as  those  of  the  plate  at  those  plate  edges  which  coincide  with 
the  beam  ends. 

1.  Separable  Eigenfunctions 

Sinusoidal  free  motion  of  a  plate  at  frequency  oi  may  be  described 
by  w(x,y)eif‘^t,  vjhere  w(x,y)  denotes  the  deflection  shape  and  obeys  the 
classical  equation  (refs.  21,  27). 

DV^w  -  =  0  ,  (D.l) 

where  D  denotes  the  flexural  rigidity  of  the  plate,  p,  its  mass  per 
unit  area.  For  a  given  set  of  plate  boundary  conditions  one  may  find 
a  doubly  infinite  set  of  solutions  w  =  iJ'jjj^(x,y)  that.  In  conjunction 
with  CD  =  05  ,  satisfy  Eq„  (D.l)  and  the  prescribed  boundary  conditions. 

Forced  motion  of  the  plate  may  be  described  by  the  addition  of  a 
term  p(x,y)  on  the  right-hand  side  of  Eq.  (D.l),  The  plate  response 
w(x, y)  and  forcing  pressure  distributions  may  be  expressed  in  terms 
of  eigenfunction  series* 


*A11  double  summations  extend  over  m,  n,  =1,2, 
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F(y) 


w(x,y)  V  Y 

v.'hose  coefficients  are  related  as  (ref.  28) 


"’mn  '  (“'2) 


w. 


mn 


(D.3) 


if  damping  of  the  plate  is  also  taken  into  account  by  replacing  the 
real  parameter  D  of  Eq.  (D.l)  by  the  complex  parameter  D*  =  D(l  +  iQ ) 
6  denotes  the  loss  factor  of  the  plate. 


A  useful  special  case  occurs  if  the  eigenfunctions  can  be 
expressed  as  products  of  separate  functions  of  the  two  plate  coordi¬ 
nates: 


For  the  usual  ideal  boundary  conditions  (clamped,  free,  or  supported) 
the  ii  are  orthogonal;  hence  the  factor  functions  are  also  ortho™ 
gonal,  and  one  may  write 


i 

i 


fni(x)  H 


Ub^  for  m  m' 
m 


0  for  m  ^  m’ 


gj^(y)  g^,  (y)  dy 


aP^  for  n  =  n' 


to  for  n  n’ 


(D.5) 


where  and  P  are  constants,  and  the  indicated  integrations  are 
carrled’"out  over  the  plate  edge  lengths. 


2.  Beam-Flate  Interaction;  Excitation  at  Beam 

The  force  distribution  Q(y)  acting  between  the  beam  and  the 
plate  may  be  Interpreted  as  a  pressure  distribution  p(x, y)  = 

Q(y)  A(x),  where  A  denotes  the  Dirac  delta  function.  The  plate- 
deflection  due  to  Q(y)  may  then  be  obtained  by  application  of  the 
foregoing  result.  In  particular,  the  plate  deflection  at  the  beam 
(at  X  =  O),  which  is  also  the  beam  deflection  and  is  denoted  here 
by  z(y)j  is  found  to  be  given  by 


z(y)  =  £ 

n 


'S 


Z  =  -r  (K  -iA  ) 
n  ixL  ^  n  n' 


(D.6) 
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where  (analogously  to  ref.  29), 


_  f^(o)  h/ai^„ 

V*  m '  ^  '  mn 

^  =  L  ^ 

m 


p  p  ^ 

m  h  -fS 


^n  =  E 


f^(o) 

m  mn  ^ 


T 


m 


P  2 
m  n  4€ 


h  K  1  -  (a)/u>^)' 


and  is  obtained  from  the  expansion  Q{y)  =  ^  Q^gj^(y). 


n 


The  beam  deflection  z,  however,  must  satisfy  the  beam’ s  equa.- 
tlon  of  motion  (refs.  21,  27) 


=  P(y)  -  Q(y)  (D.8) 

dy 

where  B*  =  B(l  +  ip)  denotes  the  beam's  (complex)  flexural  rigidity, 
M  Its  mass  per  unit  length,  and  where  F(y)  denotes  an  externally 
applied  force  distribution  as  shown  in  Pig.  14.  If  this  external 
force  is  also  expanded  as  F(y)  =  2  Pj^gj^(y),  one  may  obtain  the  fol- 

n 

lowing  relation  between  modal  force  and  modal  beam  deflection  Z^: 


n 


M 


(i+ip)  -  aj2 


K  + 

■i-  ijlL  ^ 

K  + 

n  n 


(D.9) 


.th 


Here  cs  denotes  the  beam's  n  natural  frequency,  which  may  be  shown 
to  obey 


2  „  B  ^ 

“n  “  ^n  M  '  ” 


^  «n  dy 


n 

* 


(D.IO) 


where  is  a  constant  [as  a  consequence  of  Eq.  (D.4)] 


In  absence  of  damping  (p  =  6  =  0),  infinite  amplitudes  Z  result 
for  finite  driving  force  P^  (l.e.,  resonance  occurs j  whenever^ 


Re 


^n/^n} 


0  in  Eq.  (D,9).  Damping  of  practically  feasible  magnitudes 
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generally  changes  the  resonance  frequencies  only  slightly,  so  that 
for  all  practical  purposes  one  may  still  assume  that  the  real  part 
of  Eq,  (D.Q)  vanishes  at  resonance.  Then,  In  the  presence  of  damp 
Ing,  at  resonance 


Res 


+ 


UIA 


n 


K^+A^ 
n  n 


(D.ll) 


Of  the  terms  on  the  right-hand  side  of  Eq.  (D.ll)  the  first  repre¬ 
sents  the  damping  effect  of  the  beam,  the  second  that  of  the  plate. 
From  Eqs.  (D,7)  one  finds  that  significant  contributions  to  A  occur 
only  for  ui  ™  u)  .  The  optimum  damping  effect  of  the  plate  on  oscil¬ 
lations  driven  at  the  beam  hence  may  be  realized  If  one  designs  the 
beam-plate  system  so  that  its  resonances  coincide  with  those  of  the 
plate,*  For  such  an  optimized  system  and 

('  ^  \ 

Mp+nL-^6  (D.12) 

For  many  (rectangular)  systems  one  may,  at  least  In  theory,  design 
the  plate  so  that  It  has  a  resonance  wherever  the  beam  has  a  reso¬ 
nance.  For  a  simply  supported  system  with  a  centrally  located  beam 

one  finds  that  m  =  o  ^  if 

n  mn 


1  + 


(D.13) 


In  order  for  a  plate  resonance  to  correspond  to  each  beam  resonance 
the  foregoing  equation  must  be  satisfied  for  each  value  of  n,  n  = 
1#2,3#'**.  One  thus  needs  merely  to  select  some  design  value  a  for 
the  ratio  m/n,  then  choose  the  beam  and  plate  dimensions  according 
to  Eq.  (D.13).  (If  Eq.  (D.13)  is  to  hold  for  all  n,  a  must  be  an 
Integer. } 


♦Heckl  (ref.  8)  obtained  this  conclusion  by  considering  the  plate 
as  an  array  of  rods  and  demonstrated  its  validity  experimentally. 
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3.  Excitation  on  Plate 


In  the  presence  of  an  externally  applied  pressure  r(x,y)e^“^ 
acting  on  the  plate j  one  may  write  the  modal  loading  P  of  Eq.  (D.2) 
as 


mn 


R 


mn 


+  ~m — 

m 


(D.14) 


where  the  first  term  denotes  coefficients  of  r(x,y)  and  the  second 
denotes  the  contribution  of  the  beam-plate  interaction  force  Q(y)  = 
n  From  Eq.  (D.8)  with  F(y)  =  0  one  may  establish  a  rela¬ 

tion  between  Q  and  the  modal  beam  deflection  Z  ;  the  latter  also 
obeys  ^  " 


n 


Y, 


m 


f  (0) 

mn  m'  ^ 


(D.15) 


in  view  of  the  definition  of  modal  deflection  Wj^j^  of  Eq.  (D.2). 
Combination  of  these  results  with  Eq.  (D.3)  then  permits  one  to 
arrive  at  the  following  relation  between  modal  pressure  R  and 
modal  deflection 

mn 

^  [(l+iP)a.®-®®]  X 

mn 

s 

,  M 

Equation  (D.l6)  represents  an  infinite  set  of  simultaneous 
equations.  Solution  is  difficult  since  the  equations  are  coupled, 
as  evident  by  the  appearance  of  all  W  terms  on  the  right-hand  side. 
Fortunately,  however,  one  need  not  attempt  a  solution  to  discuss 
resonant  behavior  of  the  system. 


"sn 

mn  m'  * 
(D.16) 


At  resonance  of  the  m,  n  mode  the  modal  amplitude  may  be 
expected  to  exceed  all  others  considerably.  System  resonances 
may  be  expected  to  occur  for  Re/R__/W_„l  =  0,  or  at  frequencies 

fimn  ^  J 


CJUL 


a 


mn 


m  n 


mn 


1  +  J, 


m 


(D.17) 
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The  resonant  modal  displacement  then  obeys 


R 


mn 


IxW, 


mn 


+  P'^m  “n 


Res 


(D.18) 


The  first  term  represents  tlie  plate  contribution,  the  second  that  of 
the  beam;  J  may  be  seen  to  describe  essentially  the  coupling  of  the 
n-mode  of  tne  beam  to  the  m, n-mode  of  the  plate . 

■4.  Conclusions 


Beam  resonant  responses  may  be  reduced  by  damping  the  attached 
plates.  Optimum  effective  beam  damping  is  obtained  when  the  plate 
resonances  are  made  to  coincide  with  the  beam  resonances.  For  some 
conf igui’ations  this  matching  may  be  accomplished  for  all  beam  reso¬ 
nances  by  suitable  adjustment  of  the  beam  and  plate  geometries. 


On  the  other  hand,  beam  damping  has  generally  less  effect  on 
the  resonant  responses  of  plates.  Only  a  stiff,  massive,  and  highly 
damped  beam  can  have  rtuch  influence  on  the  resonant  responses  of  a 
plate.  But  even  such  beams  have  little  effect  on  those  modes  for 
which  the  beam  Is  located  near  a  nodal  line  of  the  plate,  or  for 
which  plate  flexure  perpendicular  to  the  beam  predominates. 


E.  Damping  of  Aircraft  Structural  Joints 

As  discussed  in  the  first  section  of  this  report,  the  damping 
at  Joints  may  be  important  in  controlling  the  vibrations  of  a  struc¬ 
ture.  Increases  in  the  damping  at  Joints  might  be  obtainable  with 
relatively  minor  structural  modifications,  which  might  thus  provide 
a  simple  and  economical  means  for  reducing  vibratory  response  of  the 
structure.  This  section  summarizes  the  results  of  an  exploratory 
study  aimed  at  gaining  some  understanding  of  the  mechanisms  responsi¬ 
ble  for  the  damping  in  riveted  and  similar  Joints;  such  understanding 
is  clearly  necessary  for  rational  optimization  of  this  damping. 


Some  studies  of  similar  problems  have  pointed  toward  the  impor¬ 
tance  of  Coulomb  or  dry  friction  damping  in  some  structural  Joints. 
For  example  Klumpp  and  Goodman  (ref.  30T  developed  analyses  for 
press-fit  Joints,  and  Plan  (ref.  31)  investigated  built-up  beams. 
These  studies  showed  good  agreement  between  experimental  measurements 
and  theoretical  predictions  based  on  Coulomb  damping:  strong  depen¬ 
dences  of  damping  on  amplitude  (i.e.,  nonllnearlties)  were  found. 
Someivhat  different  results  were  reported  by  Mead  (ref.  32).  In  ex¬ 


periments  performed  on  more  nearly  practical  riveted  Joints  he 
observed  no  amplitude  dependence  of  damping,  except  at  relatively 
high  amplitudes  (which  are  of  lesser  Interest  here,  since  one 
desires  to  avoid  severe  vibrations  in  the  first  place). 


kk 
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As  Mead  has  pointed  cut^  the  damping  of  riveted  joints  might  be 
affected  by  many  factors^  Including: 

1.  Amplitude  and  mode  shape  of  the  vibration. 

2.  Normal  pressure  between  Joined  structures  and  Its  distribution 
about  rivets. 

3.  Shear  stiffness  of  rivet  in  relation  to  local  stiffness  of 
Joined  structures.  End  fixity  of  rivets  associated  with  rivet 
head  geometry^  radial  pressure  between  rivet  and  hole;  l.e., 
assembly  methods,  might  make  significant  differences. 

4.  Surface  condlt.ions  at  the  interface.  Including  roi;:ghness,  presence 
of  lubricants,  dust,  local  oxidation,  wear  (fretting). 

In  order  to  assess  the  Importance  and  trends  associated  with  some 
of  these  parameters,  a  number  of  damping  measurements  were  undertaken 
using  essentially  the  experimental  setup  described  by  Kerwln  (refs. 

33t  34).  The  measurements  consisted  of  driving  a  test  bar  sinusoidally, 
removing  the  excitation,  and  observing,  the  decay  of  the  vibration.  Prom 
an  oscilloscope  display  of  the  logarithm  of  amplitude  versus  time  the 
linearity  of  the  decay  could  be  assessed  directly,  and  the  loss  factor 
could  be  obtained  from  the  slope  of  the  decay  curve  by  means  of  a  simple 
calculation  (ref.  34). 

1.  Experimental  Results 

A  preliminary  experiment  was  performed  on  two  1/8"  thick  24ST 
aluminum  plates,  2"  wide  and  15.5"  lor^g#  Joined  by  8  bolts  (arranged 
4"  apart  in  two  parallel  rows,  with  l"  between  rows).  The  decay 
curves  were  observed  to  be  straight,  l.e.,  damping  was  found  to  be 
linear,  over  the  entire  range  of  amplitudes  (0.02  to  50  g)  the 
experiment.  Similar  linearity  was  also  obsei’ved  In  virtually  all  of 
the  subsequently  reported  experiments,  but  In  most  of  these  no  attempt 
was  made  to  obtain  a  wide  range  of  amplitudes. 

A  series  of  damping  measurements  were  carried  out  on  test  pieces 
cut  from  the  interior  of  the  forward  structure  of  an  F-105  aircraft. 
These  test  samples,  whose  cross-sections  are  sketched  In  Fig.  15,  con¬ 
sisted  of  27  Inch  long  beams  with  sections  of  skin  attached  by  rivets 
or  spot  welds.  Test  results  for  these  samples  appear  In  Fig.  15,  to¬ 
gether  with  results  pertaining  to  the  riveted  sample  piece  sifter  the 
rivets  had  been  drilled  out  and  replaced  by  bolts  and  nuts.  From 
these  results  It  appears  plausible  that  the  different  Joining  methods 
considered  result  in  essentially  the  same  damping  If  the  geometries 
of  the  Joints  are  roughly  similar.  Since  Fig,  Ip  pertains  to  methods 
that  produce  tight  Joints,  l.e.,  that  produce  good  connections  in  the 
vicinity  of  the  connectors,  the  foregoing  statersent  might  be  extended 
to  include  any  tight  Joint.  It  thus  appears  that  the  damping  of  such 
tight  Joints  is  not  due  to  the  connectors  themselves,  but  perhaps  due 
to  relative  motion  at  the  Interface  at  some  distance  from  the  connec¬ 
tors. 
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In  view  of  Fig,  15^  one  may  also  expect  to  be  able  to  obtain 
results  that  have  some  meaning  in  relation  to  riveted  Joints  fi’om 
experiments  conducted  on  bolted  Joints,  which  are  much  easier  to 
adjust  to  desired  conditions.  Figure  16  marshals  additional  data 
In  support  of  the  similarity  of  the  damping  characteristics  of 
riveted  and  bolted  Joints.  (No  torque  wrench  was  available  for 
the  small,  4-40,  screws  that  were  used.  The  "tight”  condition 
corresponds  to  the  tightest  that  could  be  obtained  by  hand  with 
a  screwdriver;  the  "loose"  was  obtained  by  backing  the  screws  off 
one  quarter  turn  from  the  tight  condition.) 

Effects  of  rivet  spacing  may  be  seen  in  Fig.  17,  which  sum¬ 
marizes  the  damping  results  obtained  for  the  aforementioned  riveted 
aircraft  Joint  sample  before  and  after  drilling  out  some  of  the 
rivets.  Again,  since  damping  Is  certainly  not  proportional  to  the 
number  of  rivets,  it  cannot  be  due  to  the  rivets  per  se  or  to  small 
areas  surrounding  the  i*ivets.  It  may  be  more  reasonable  to  asso¬ 
ciate  damping  with  relative  motion  at  interface  areas  located  out¬ 
side  of  "clamped"  areas  near  connectors.  From  the  near  correspon¬ 
dence  of  the  curve  for  38  to  that  for  I8  rivets  one  might  estimate 
that  in  the  original  sample  the  clamped  areas  associated  with  ad¬ 
jacent  rivets  probably  overlapped. 

That  damping  can  not  be  associated  merely  with  the  structure 
In  the  Immediate  neighborhood  of  the  connectors  is  also  evident 
from  Fig.  18.  The  curves  of  this  figure  pertain  to  the  spot-welded 
sample  structure,  with  different  widths  of  the  skin  portions  (cor¬ 
responding  to  the  3”  dimension  of  Fig.  I5).  Since  the  Identical 
connections  are  involved  in  all  the  cases  shown  here;  they  hence 
cannot  be  responsible  for  the  changes  observed,  and  one  is  once 
again  led  to  ascribe  the  damping  to  those  portions  of  the  struc¬ 
ture  which  are  not  near  the  spot  welds. 

Figures  19  and  20  Illustrate  the  effects  of  Interface  lubri¬ 
cation.  The  addition  of  lubricant  was  found  to  result  In  Increased 
damping,  with  less  viscous  lubricants  producing  higher  damping.* 
Comparison  of  these  two  figures  also  indicates  that  reduction  of 
interface  pressure  (loosening  of  the  screws)  results  in  increased 
damping  for  all  of  the  interface  conditions  tested.  Pressure  reduc¬ 
tion  and  improved  lubrication  thus  produce  similar  results,  probably 


*It  should  be  noted  that  some  lubricant  was  found  In  the  Joint  as 
cut  from  the  aircraft;  the  corresponding  damping  curve  is  labeled 
"original  condition".  Chemical  degreasing  resulted  in  reduced 
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coating  (curve  labeled  cleaned  to  metal)  lowered  the  damping  even 
more . 
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FIG.  16  COMPARISON  OF  RIVETS  AND  VARIOUSLY  TIGHTENED  SCREWS 
(SAMPLE  NO. 2) 
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F1G.I8  EFFECTS  OF  WIDTH  OF  ATTACHED  SKIN 
PORTION  (SAMPLE  NO.  I  ) 
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F1G.I9  EFFECTS  OF  INTERFACE  LUBRICATION 
(SAMPLE  N0.2  WITH  18  SCREWS,  TIGHT) 
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FIG.20  EFFECTS  OF  INTERFACE  LUBRICATION 
(SAMPLE  N0.2  WITH  18  SCREWS,  LOOSE) 
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since  both  serve  to  reduce  Interface  friction.  Reduced  friction 
results  in  increased  relative  motion,  which  evidently  more  than 
compensates  for  the  effect  of  the  reduced  friction  as  far  as  energy 
dissipation  is  concerned.  (Note  that  the  cyclic  energy  dissipated 
in  a  classical  linear  dashpot  is  proportional  to  the  square  of  the 
amplitude,  but  only  to  the  first  power  of  the  viscous  damping  co¬ 
efficient.  ) 

Figure  21  shows  the  results  that  some  simple  modifications  had 
on  the  damping  of  the  riveted  P-I05  sample  structwe.  With  the  toe 
of  the  original  equilateral  L-shaped  cross-section  bent  down  a 
general  Increase  in  damping  was  observed,  which  may  at  least  in  part 
be  attributed  to  the  extremely  good  contact  fostered  by  such  a  cross- 
section  at  the  heel  and  to«j  of  the  L.  However,  a  further  general 
increase  in  damping  was  obtained  (except  below  630  cps)  when  the  toe 
of  the  L  was  cut  off  altogether.  This  increase  might  be  due  to  im¬ 
proved  contact  in  the  area  that  previously  was  between  the  heel  and 
toej  such  Improvement  would  be  most  noticeable  at  the  higher  fre¬ 
quencies,  where  the  motion  is  less  restricted  by  the  "rivets".  The 
peak  observed  at  about  500  cps  with  the  bent-down  toe  is  probably 
due  to  a  mode  shape  for  which  a  considerable  amount  of  relative 
motion  occurs  at  the  toe. 

Figures  22  and  23  summarize  results  of  experiments  aimed  at 
evaluating  the  effects  of  beam-and-skln  contact  area  and  of  beam 
stiffness  separately.  The  measurements  indicated  in  these  figures 
were  performed  on  specially  designed  aluminum  channel  beams;  those 
of  Fig.  22  have  nearly  the  same  flexural  rigidity  but  differ  con¬ 
siderably  in  contact  area,  whereas  those  of  Pig.  23  have  eq;;ial  con¬ 
tact  areas  but  different  stiffnesses. 

Increased  contact  area  was  found  to  result  In  increased  damping 
and  changes  in  beam  stiffness  were  found  to  affect  damping  relatively 
little  in  general.  The  picture  is  by  no  means  clear,  however.  It  is 
difficult  to  say,  for  example,  if  the  higher  damping  noted  with  in¬ 
creased  contact  area  is  due  to  this  area  Increase  per  se  or  to  changes 
in  mode  shape.  The  high  peaks  shown  in  Fig.  23  between  5OO  and  800  cps 
also  present  somewhat  of  an  anomaly.  Calculations  show  that  this  fre¬ 
quency  range  corresponds  to  that  where  one  half  wavelength  of  flexural 
motion  of  the  akin  is  equal  to  the  width  of  the  skin  attached  to  the 
test  beams,  so  that  these  peaks  might  be  due  to  bending  of  the  com¬ 
posite  about  the  rivet  line  (in  contrast  to  the  rest  of  the  curve 
which  probably  is  primarily  associated  with  bending  of  the  rivet  line), 

A  final  series  of  tests  were  performed  to  determine  the  effects  of 
adhesives.  The  results  of  these  tests,  as  shown  in  Pig.  24,  indicate 
that  a  dissipative  adhesive  like  MMM  ^l66  can  result  in  a  considerable 
damping  Increase  over  a  wide  frequency  range,  whereas  a  relatively 
rigid  loss-less  adhesive  (epoxy)  not  only  adds  little  damping  but  even 
reduces  the  damping  of  riveted  Joints  under  some  circumstances. 
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FIG.2i  EFFECTS  OF  CONTACT  AREA  MOD  IF  !G  AT  IONS  OF  SAMPLE  N0.2. 
(CLEAN  INTERFACE,  18  SCREV/S,  TIGHT  ) 


FtG.22  EFFECT  OF  CONTACT  AREA  ON  DAMPING  OF 
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EFFECT  OF  BEAM  STIFFNESS  ON  DAMPING 
OF  BOLTED  JOINT. 
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•  Conclusions 


The  studies  outlined  here,  althoijgh  by  no  means  conclusive  or 
as  comprehensive  as  one  would  like,  permit  one  to  arrive  at  least 
at  some  tentative  conclus  ions  * 

These  are: 

1.  The  damping  of  aircraft  joints  is  essentially  linear,  except 
perhaps  at  high  amplitudes. 

2.  The  damping  of  tight  Joints  seems  not  to  be  primarily  associated 
with  the  connectors  or  the  structure  in  the  immediate  vicinity 
of  the  connectors. 

a.  Different  connectors  (spot  weld,  rivets,  bolts)  produce 
similar  damping  characteristics. 

b.  Studies  on  bolted  joints  are  capable  of  yielding  results 
that  have  some  meaning  in  relation  to  riveted  joints. 

3.  Looser  joints  produce  higher  loss  factors;  reduced  interface 
pressure  (looser  or  fewer  rivets)  or  better  lubrication  results 
in  higher  damping. 

4.  Increased  beam-to-skin  contact  area  produces  increased  damping. 

5.  Beam  stiffness  per  se  has  relatively  little  effect  on  damping. 

6.  Rubbery  dissipative  adhesives  introduced  in  the  joint  can 
increase  damping  drastically,  but  rigid  lops-lesa  adhesives 
can  reduce  damping. 

Although  the  results  presented  here  apply  to  a  s‘*-ructural  Joint 
configuration  which  appears  to  be  fairly  typical  of  current  practice, 
one  should  use  extreme  care  in  extrapolating  these  results  to  other 
structures.  Available  data  are  still  too  limited  and  understanding 
of  the  pertinent  phenomena  is  still  to  incomplete  to  permit  one  to 
draw  valid  conclusions  about  the  damping  of  Joints  of  structures  in 
general . 

It  should  also  be  noted  that  the  various  high  peaks  observed  in 
the  presented  damping  curves  are  probably  associated  with  modal  shapes 
that  favor  high  energy  dissipation  In  some  manner.  Since  such  shapes 
are  probably  altered  or  suppressed  if  the  Joint  is  part  of  a  more 
complicated  structure,  such  peaks  are  most  likely  not  observed  in  such 
structures.  In  judging  the  damping  effectiveness  of  a  joint  which  is 
to  be  incorporated  in  a  more  complicated  structure  from  damping 
measurements  performed  on  the  joint  by  itself  one  should  therefore 
probably  ignore  the  peaks  and  concentrate  on  the  overall  level  of 
damping. 
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F.  R€^lation  of  Plate  Absorption  Coefficients 
to  Loss  Factoids  of  Boundar'y  Structures 

A  considerable  amount  of  information  is  available  concerning 
the  loss  factors  of  beam-like  structures,  whereas  very  little  is 
known  about  absorption  coefficients  associated  with  such  structures. 
This  section  tlierefore  is  concerned  with  some  preliminary  analytical 
steps  toward  obtaining  interrelations  between  these  two  quantities. 

1.  general  Expressions 

Prom  the  first  section  of  this  report  it  is  evident  that  the 
absorption  coefficient  7  of  a  portion  of  a  plate  boundary  may  be 
defined  by 


“  ""^diSB^/^^inc^  ^ 

where  the  (length-wise)  average  power  flow  per  unit 

length  measured  along  the  crest  of  a  plate  wave  incident  on  the 
boundary.  Pdiss  <ienotes  the  (length-wise)  average  power  dissipated 
at  the  boundary,  again  referred  to  unit  length  measured  along  the 
incident  wave  crest.  The  brackets  <  >  indicate  averages  taken  over 
the  angle  of  incidence  0,  so  that  7  denotes  the  fraction  of  the 
Incident  povfer  that  is  dissipated. 

If  the  wave  field  on  the  plate  Is  reverberant  (l.e,.  Independent 
of  0)  then  (ref.  4) 

<Plj^c>  ^  ^Inc  =  5Cp  D  k|  0)  ,  (F.2) 

where  Xp  denotes  the  displacement  amplitude  of  the  plate  flexural 
wave,  kp  the  wave  number,  tu  the  circular  frequency,  and  D  the 
flexural  rigidity  of  the  plate. 

The  loss  factor  q  of  a  structure  is  defined  by  (ref.  2) 

where  Ujj[ss  denotes  the  total  energy  dissipated  by  the  structure 
per  cycle  and  Uq  the  "energy  of  vibration"  of  the  structure.  (For 
lightly  damped  structures  Uq  is  essentially  equal  to  the  time-wise 
maximum  of  the  total  strain  energy).  In  general  n  for  a  given 
structure  depends  on  amplitude,  mode  shape,  and  frequency.  No 
essential  difficulty  results  if  the  energy  terms  of  Eq.  (F.3)  are 
intei-pi-eted  as  length-v,'ise  average  values  rather  than  totals;  hence 
such  averages  will  be  assumed  henceforth.  If  q  is  now  assumed  to 
be  Independent  of  length-wise  position  in  the  structure,  then 
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’’diss  =  ''dlss  -  U^r,<» 


(F.4) 


and  from  Eq.  (P.l) 

Y  =  <f\^Q>Ap  °  (F-5) 

2.  Boundary  Structures  with  Lobs  Factors 

Independent  of  Deflection  Shape _ 

For  beam-like  structures  whose  loss  factors  do  not  depend  on 
the  shapes  into  which  they  are  deformed  during  a  vibration  (e.g., 
homogeneous  viscoelastic  bars  or  elastic  bars  with  free  visco¬ 
elastic  layers),  the  energy  losses  obtained  when  such  structures 
are  attached  to  plates  will  be  independent  of  the  angle  of  -inci¬ 
dence  of  plate  waves,  and  <riU^>  =  'n<U^>, 

The  average  strain  energy  in  a  beam  of  flexural  stiffness  El 
(where  one  might  consider  an  ^'effective  beam"  as  including  a  por¬ 
tion  of  the  plate)  may  be  found  to  be  given  by  (ref.  4) 

<V  “  ^ 

If  the  beam  Is  assumed  to  deflect  sinusoidally  (in  space)  with  the 
proper  spatial  periodicity  to  conform  to  the  Incident  plate  waves 
(refs,  4,  35)*  From  Eq.  (F.3)  one  may  then  obtain 

7  =  ti(EI  kp/4D)-R  ;  R  s  <(XgAp)^  sin^0>  .  (P.7) 

The  parameter  R  depends  on  the  effect  of  the  beam  on  the 
plate  motion.  When  the  added  beam  is  so  light  and  soft  as  to  have 
no  appreciable  effect  on  this  motion,  then  Xg«Xp  and  R»3/B,  In 
general,  however,  more  detailed  analyses  are  required.  Some 
results  directly  applicable  to  such  analyses  are  available  (ref, 
35)j  tint  the  multitude  of  parameters  complicate  the  analysis  to 
such  an  extent  that  it  must  be  considered  beyond  the  scope  of  the 
present  study. 

3.  Boundary  Structures  Attached  with  Viscoelastic  Adhesive 

In  the  case  where  structural  beam  is  Joined  to  a  plate  by  a 
thin  layer  of  viscoelastic  glue,  virtually  all  of  the  energy  dis¬ 
sipation  occurs  in  the  viscoelastic  material.  One  may  then  utilize 
the  definition  of  the  material's  shear  loss  factor  to  write 

Udlss  =  ,  (P.8) 
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where  Ug  denotes  the  tlme-vjlse  maximum  of  the  (length-wise  average^ 
strain  energy  stored  in  the  viscoelastic  layer.  Fi’om  Eqs.  (F.l  -  4) 
and  the  foregoing  one  may  then  obtain  the  alternate  expression 

7  =  KUg>Ap  D  (F.9) 

Prom  an  analysis  similar  to  Kerwin'a  (ref.  33)  one  may  show 
that  for  a  laminate  consisting  of  a  thin  viscoelastic  layer  between 
two  elastic  structures  (i.e.^  a  beam  and  plate)  where  these  are  much 
stiff er  (in  extension)  than  the  viscoelastic  layer^ 


“b  =  TO  ‘S  Xg 


sln^& 


y  +  sln^0 


(F.IO) 


if  one  also  assumes  that  the  wave  in  the  laminate  is  sinusoidal  in 
space  and  has  its  wavelength  determined  by  plate  waves  striking  the 
beam  with  angle  of  incidence  0.  The  newly  Introduced  symbols  have 
the  following  meanings;* 


y 


(F.ll) 


b  =  width  1 


H  = 
G  = 
S  = 


thickness 

real  part  of  complex  shear 

distance  from  neutral  axis 
that  of  added  beam 


modulus 


of  viscoelastic 
layer 


of  plate  to 


Kg  =  EgAs  =  extenslonal  stiffness  of  unit 
length  of  added  beam 

?■  =  Epbph  =  extenslonal  stiffness  of  unit 

length  of  strip  of  plate  affected  by  beam 

Eg,Ep  =  Yoiang's  modulus  of  beam,  plate,  material 

rtg  =  beam  cross-sectional  area 

h  =  plate  thickness 

bp  =  effective  width  of  plate  strip 


*The  parameter  y  is  closely  related  to  the  generalized  shear  parameter 
of  Section  III  of  Ref,  2, 
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h ,  Damping:  Tapes 


If  X-  is  indepcpdent  of  the  angle  of  Incidence  6,  then  one  may 
readily  average  Ug  of  Eq,  (P.IO)  over  Q.  Introduction  of  this  result 
into  Eq.  (F.9)  then  permits  one  to  arrive  at 


(F.12a.) 


where 


! 


IL 

2 


4y  + 


(P.12b) 


Damping  tapes  represent  the  simplest  boundary  structures  for 
which  these  relations  apply.  They  are  studied  here  not  only  because 
of  their  possible  practical  importance,  but  also  because  such  a 
study  might  provide  an  indication  of  the  validity  of  the  foregoing 
analysis . 


Figure  25  shows  absorption  coefficients  y  measured  with  2-lnch 
wide  strips  of  damping  tape  applied  to  an  Irregularly  shaped  aluminum 
plate,  as  well  as  curves  obtained  from  calculations  based  on  Eqs, 
(F.I2).  It  is  evident  that  theory  and  experiment  agree  quite  well  as 
long  as  the  tape  width  is  less  than  the  half  wavelength  of  a  plate 
flexural  wave,  l.e.,  for  X<4  Inches.  The  derivation  of  Eqs.  (P.12) 
assumes  that  all  damping  is  due  to  the  length-wise  deformation  of  the 
viscoelastic  layer.  This  assumption  Is  no  longer  adequate  for  higher 
frequencies  where  the  tape  width  may  extend  over  several  plate  wave¬ 
lengths.  At  these  frequencies  one  would  expect  to  obtain  additional 
energy  dissipation  due  to  width-wise  shear  deformations  of  the  visco¬ 
elastic  layer,  and  Eqs.  (F,12)  underestimate  the  absorption  coeffi¬ 
cient.  This  is  also  evident  in  tne  figure. 


The  experimental  v/ork  summarized  in  Pig.  25  was  done  by  methods 
discussed  in  Section  I  and  reference  3^.  Commercial  (Scotch  Brand) 
damping  tapes  were  used.  The  thicknesses  of  the  adhesive  layers  of 
these  tapes  are  a  nominal  5  mils,  428B  tape  has  8  mil  aluminum  foil 
backing,  428C  has  a  similar  12  mil  thick  backing.  The  adhesive  prop¬ 
erties  used  in  the  calculations  were  those  published  in  reference  33- 
For  the  calculations  Xjj^Xp  was  asBujned,  since  the  tape  is  much  lighter 
and  softer  than  the  plate  and  may  be  supposed  to  affect  the  plate 
motion  only  slightly. 
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•/is"  thick  aluminum  plate 


SECTION  III 


SUr-lMARy  AND  DISCUSSION 


It  has  been  shown  that  the  responses  of  complex  structures 
are  strongly  Influenced  by  energy  absorption  at  joints  between 
substructures  (measured  in  terms  of  absorption  coefficients)  and 
by  the  amount  of  coupling  between  an  exciting  sound  field  and 
the  structure  (described  in  terms  of  radiation  resistance )j  in 
addition  to  the  mox*e  conventional  mass  stiffness  and  damping 
parameters  of  vibration  analysis. 

Absorption  coefficients  were  measured  for  a  few  structural 
Joints,  primarily  to  validate  some  of  the  theoretical  work.  De¬ 
tailed  studies  of  absorptive  etructui’es  and  of  some  of  the  energy 
dissipation  meclianlsms  associated  with  them  remain  to  be  under¬ 
taken,  particularly  with  a  view  toward  developing  configurations 
and  design  data  leading  to  improved  absorption  performance. 


The  dsunplng  due  to  riveted  and  similar  Joints,  which  now 
appears  to  be  much  more  significant  than  had  previously  been  sup¬ 
posed,  should  merit  considerable  further  investigation,  Some  of 
the  features  of  Joint  damping  bro\jght  to  light  by  the  exploratory 
study  reported  here  should  be  explored  further]  their  exploitation 
seems  likely  to  lead  to  damping  increases  secured  at  little  expense 
and  weight  penalty. 

An  Introductory  analytical  study  attempted  to  establish  some 
relations  between  loss  factors  of  boundary  structures  and  the  ab¬ 
sorption  coefficients  realized  when  these  structures  are  attached 
to  plates.  Since  much  loss  factor  data  are  available,  whereas 
relatively  little  is  known  about  absorption  coefficients,  continu¬ 
ation  of  these  studies  may  lead  to  results  of  considerable  utility. 


Reductions  in  radiation  resistance  result  in  reduced  vibra¬ 
tions  of  structures  exposed  to  acoustic  fields,  hence  control  of 
these  resistances  constitutes  an  Important  means  of  vibration  con¬ 
trol.  The  first  steps  toward  providing  an  understanding  of  the 
mechanisms  responsible  for  radiation  resistance  were  described  in 
this  report]  extensions  of  these  steps  are  being  studied  (ref.  17 )> 
and  undoubtedly  much  additional  investigation  will  be  required 
before  generally  useful  prediction  and  design  procedures  applicable 
to  complex  structures  are  developed.  The  ability  of  special  beams 
(designed  per  Reference  19)  to  reduce  radiation  resistance  has  been 


demonstrated,  h 

A  A"'!  -iiTi 


owever  only 


by  means  of  a  laboratory  experiment. 

_ _ i _ _  ‘  ~  z  r  -r  _.._u 
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structures  is  warranted  and  feasible  in  practice. 


ASD  TDR  62-237 


64 


Detailed  analytical  studies  of  vibration  absorbers  utilizing 
viscoelastic  spring  elements  showed  these  absorbers  to  be  rela*' 
tively  Impractical.  The  same  conclusion  was  reached  concerning 
the  utility  of  absorbera  using  distributed  systems. 

Viscoelastic  leaf  springs  were  found  to  hold  some  promise 
as  vibration  isolation  systems  required  to  perform  well  over  very 
wide  frequency  ranges.  Their  main  advantages  over  more  conventional 
springs  lie  in  that  they  exhibit  fewer  transmlsslblllty  peaks  and 
that  they  can  be  designed  with  frequency- variable  damping.  The 
latter  is  Important  in  reducing  the  magnitude  of  the  transmlssibility 
peaks  without  increasing  significantly  the  transmlssibility  at  low 
frnquoncies.  Additional  study  la  needed  to  reduce  these  concepts 
to  practical  hard-ware, 

A  study  of  the  modal  behavior  of  beam-plate  systems  showed 
that  when  such  systems  are  excited  at  a  beam  the  plates  adjacent 
to  the  beam  may  contribute  significantly  to  I'eduction  of  the  reso¬ 
nant  beam  responses.  High  effective  beam  damping  is  obtained  when 
the  plates  are  highly  damped  and  when  they  ars  designed  so  that 
they  have  resonances  whenever  the  beam  has  resonances. 

On  the  other  hand,  beam-plate  systems  excited  by  distributed 
loading  on  the  panels  tend  to  derive  only  minor  resonant  response 
reduction  from  dsunplng  of  the  beams.  This  is  so  because  the  beams 
tend  to  couple  poorly  to  the  plates.  (Plates  have  many  more  reso¬ 
nances  than  beams,  hence  one  can  not  match  beam  resonances  to  each 
of  the  panel  resonances.)  In  general,  therefore,  one  may  hope  to 
obtain  good  control  of  beam-plate  system  resonant  responses,  no 
matter  how  the  system,  is  driven,  by  damping  only  the  plates  and 
designing  these  so  that  the  resonances  coincide  with  beam  reso¬ 
nances. 
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APPENDIX  I 


MODAL  DENSITIES 


A.  Plates 


The  natural  frequencies  co  of  simply  supported  rectangular 
plates  are  given  by 


^TTC  j 


(I.l) 


where  n  are  positive  integers^  a,  b  are  the  edge  lengths^  and 
c  =  ^0)  K  is  the  wave  velocity,* 

If  (m/a)  and  (n/b)  v/ere  continuous  variables,  Eq.  (I.l)  would 
represent  a  circle  of  radius  (o/ttc)  in  the  plane  of  these  variables, 
as  shown  in  Pig.  26.  However^  the  natural  frequencies  correspond 
only  to  Integral  values  of  m  and  n,  hence  to  the  "lattice  points" 
indicated  in  the  figure.  The  number  of  natural  frequencies  that 
fall  below  a  specified  value  of  cd  thus  corresponds  to  the  number 
of  lattice  points  that  fall  within  the  quarter  circle  of  radius 
oi/irc. 


Actual  counting  of  these  points  usually  is  impractical.  It  is 
more  convenient  to  assign  an  area^i^^j  =  to  each  lattice  point 

(as  shown  in  Fig.  26)  and  to  determine  how  many  of  these  lattice 
point  areas  make  up  the  area  of  the  quarter  circle.  Of  course,  one 
must  make  due  allowance  for  the  strips  along  the  coordinate  axes 
that  are  not  associated  with  lattice  points.  One  thus  finds  that 
for  a  simply  supported  plate  the  number  of  modes  whose  frequencies 
CD  are  less  than  co  are  given  by 


N 


(cD  <<d) 
'  mn^  ' 


(1.2) 


where  P  =  2(a+b)  denotes  the  plate  perimeter.  The  modal  density, 
obtained  by  differentiation  of  (1.2)  with  respect  to  oo,  is 


ng(cD) 


1  P 


(1.3) 


*K  is  the  radius  of  gyration  of  the  cross-section,  equal  to  h/*Vi2, 
where  h  represents  the  plate  thickness,  ^  fE/p  is  the  longitu¬ 
dinal  wave  velocity  in  the  plate  material.*^  ^ 
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For  other  than  simply  supported  boundaries  the  mode  arrange¬ 
ment  is  not  quite  so  easily  computed,  but  approximations  are 
possible.  Noting  that  the  natvn’ai  frequencies  of  a  clamped-clamped 

Ui  , 

bar  of  length  b  are  approximately  given  by  ••  2,  ••  •, 

one  may  visual  i?;e  the  lattice  corresponding  to  a  plate  clamped  at 

y=0,b  to  be  shifted  upward  by  an  amount  v/.lth  respect  to  the 
lattice  of  Fig,  26.  Similarly,  clamping  at  x=0, a  v/ould  shift  the 

lattice  to  the  right  by  an  amount  ,  If  one  now  accounts  for  the 

areas  within  the  quarter  circle  that  are  not  associated  with  lattice 
points,  one  may  obtain  the  number  of  modes,  as  before.  For  a  fully 
clamped  plate  one  finds 


(1.4) 


Plates  with  free  edges  may  be  treated  simillarly.  Since  the 
natural  frequencies  of  a  bar  with  free-free  ends  are  exactly  the 
same  as  those  of  the  same  bar  with  clamped-clamped  ends,  the  lattice 
for  a  free  plate  should  be  the  same  as  that  for  a  fully  clamped  one. 
However,  a  free  plate  may  also  deform  In  only  one  coordinate  direc¬ 
tion  at  a  time,  i.e.,  like  a  bar,  so  that  additional  lattice  points 
occur  where  the  lattice  lines  intersect  the  coordinate  axes.  For  a 
free  plate,  then. 


Note  that  at  high  frequencies  one  obtains  the  same  medal 
density  for  all  boundary  conditions,  natmely 


(1.6) 
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This  corresponds  to  a  frequency  separation  between  modes  of 


B. 


Af 


&CD 

2v 


1 


Beam-Plaite  Systems 


(1.7) 


Consider  an  agglomeration  of  panels,  frames,  and  stringers, 
as  is  common  to  aircraft  construction.  If  the  frames  and  stringers 
were  infinitely  rigid  they  would  remain  motionless,  and  the  modal 
density  would  simply  be  the  sura  of  the  modal  densities  for  all  of 
the  panels.  Actual  reinforcing  structures  possess  some  flexibility, 
so  that  adjacent  panels  may  be  expected  to  interact  and  the  modes 
of  the  total  structure  will  be  complicated  combinations  of  motions 
of  all  the  sub- structures,  but  one  may  hope  to  obtain  a  rough  es¬ 
timate  of  the  modal  densities  of  such  complicated  structures  by 
assuming  that  each  panel  adds  its  modea  Independently,  The  panels 
(which  may  be  assumed  clamped  at  all  frames  and  stringers)  thus 
contribute  the  following  number  of  modes  below  the  frequency  cs; 


(1.8) 


where  v.,  denotes  the  number  of  the  panels.  For  panels  of  the  same 
material  and  of  uniform  thickness,  Eq,  (1.8)  may  also  be  written  as 


(1.9) 

(1. 10) 


where  B  is  the  total  beam  (or  rivet- line)  length.  Note  that  each 
beam  is  counted  twice  as  perimeter  length,  once  for  each  panel  it 
bounds,  making  PogB. 


[If  one  wishes  to  refine  the  estimate  one  must  consider  inter¬ 
action  of  the  various  structural  components.  Among  others,  one  may 
wish  to  add  the  modes  due  to  motion  of  the  gross  structure.  The 
motion  of  a  multiply  beam-reinforced  plate  at  the  lower  frequencies 
is  essentially  that  of  an  crthotroplc  plate  (ref,  l6),  whose  stiff¬ 
ness  is  generally  due  primarily  to  the  beams.  If  this  equivalent 


ABD  TDK  62-237 


69 


plate  la  considered  unsupported,  one  may  Invoke  Eq.  (1.5)  l>ut  use  a 
new  equivalent  phase  velocity  Cg.  If  this  orthotropic  plate  has  an 
area  equal  to  the  total  structure  area  A,  then  it  contributes  a  num¬ 
ber  of  inodes  given  by 


N®  (oj) 


(1. 11) 


below  0). 


The  number  of  inodes  below  oi  in  the  beam-plate  system  will  be 
given  by 


ca^kf  1  ^  1  \  oiB 


(1.12) 


However,  usually  c  »c,  since  the  reinforcing  beams  in  practical 
structures  generally  are  considerably  stiffen  than  the  panels. 
Thus,  Inclusion  of  the  equivalent  plate  modes  gcnoraily  modifies 
Eq.  (1,9)  only  slightly,  and  generally  adds  little  to  the  quality 
of  the  estimate.] 


The  frequency  at  which  the  first  mode  occurs  may  be  estimated 
by  setting  Nfco)  =  1  in  Eq,  (1,8),  Solution  of  the  resulting  quad¬ 
ratic  equation  results  in  a  somewhat  cumbersome  expression,  however. 
One  may  instead  wish  to  consider  when  one  obtains  one  mode  per  panel 
on  the  average,  l.e.,  vvhen  N  =  v  The  corresponding  frequency  f 
of  the  "onset  of  inodes"  then  is  obtained  from  Eq.  (1.8)  simply  as° 


(1.13) 


For  the  simulated  aircraft  panel  of  Fig.  3  Eq.  (I. 10)  results 
In  a  modal  density  value  (at  high  frequencies)  of  0,272  modes/cps, 
whereas  experimentally  a  value  of  0.28  vias  determined.  Equation 
(1.13)  predicts  onset  of  modes  (l.e.,  the  occurrence  of  N=v  -5I5 
modes)  at  I70  cps  for  the  aircraft  panel,  but  experimentally  one 
finds  N=15  at  100  cps  -  a  discrepancy  that  may  to  some  extent  be 
due  to  experimental  difficulties,  but  most  likely  is  due  primarily 
to  failure  of  the  various  assumptions  to  hold  for  the  lower  modes 
and  frequencies. 
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APPENDIX  II 


MEASURED  LOSS  FACTORS  OF  AIRCRAb'T  SUBSTRUCTTJRES 


The  loss  factors  of  some  pieces  of  a  F-IO5  for\ffard  structure 
were  measured  during  the  course  of  some  of  the  previously  reported 
studies.  These  results  are  presented  here  In  Figures  27,  and  28. 
Inclusion  of  these  data  was  deemed  desirable  here  because  such  in¬ 
formation  seems  to  be  lacking  in  the  literature  and  because  it 
might  have  some  bearing  on  the  course  of  future  damping  investiga¬ 
tions  . 

The  data  reported  here  were  obtained  by  decay  rate  measure¬ 
ments,  again  using  essentially  the  techniques  and  Instruments 
discussed  in  reference  3^. 

The  substructures  to  which  these  data  pertain  Included  no 
specific  damping  provisions.  They  were  complicated  structures, 
however,  made  up  of  many  panels  with  numerous  reinforcing  stringers 
and  bulkheads,  and  a  multitude  of  rivets,  spotwelds,  and  (relatively 
rigid)  inserts  of  various  sorts.  The  structures  were  painted  for 
the  most  part. 

Rough  measurement  of  surface  area  and  total  rivet  line  length 
showed  that  on  the  basis  of  the  relation  near  the  top  of  p.  4a 
mean  free  path  length  of  6  Inches  or  less  would  apply  for  flex¬ 
ural  waves  on  the  surface  panels  of  these  structures.  For  these 
1/16”  thick  aluminum  panels  a  6"  wavelength  corresponds  to  about 
600  cps,  so  that  the  absorption  coefficient  concept  might  reason¬ 
ably  be  applied  to  these  structures  for  somewhat  higher  frequencies. 
Unfortunately,  detailed  pertinent  rivet  line  absorption  coefficient 
data  was  lacking,  so  that  further  validation  of  the  concept  could 
not  be  accomplished.  This  validation  must  be  deferred  to  a  later 
study. 
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F  i  G.  27 


LOSS  FACTOR  OF  FRONT  END  CONE  FROM  F-105 
NOSE  SECTION 
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FREQUENCY  (CPS) 


FIG.28  LOSS  FACTOR  OF  DOOR  PANEL  FROM  F-105 
NOSE  SECTION 
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